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ABSTRACT 


After a brief introduction on the parameter plane methods, the be- 
havior of simple systems with time delay has been initially studied by 
the root-loci method and then the state of the art of parameter plane 
methods with time delay has been presented. 

A variation of the solution in the state of the art presentation 
has been derived and an appropriate digital computer program has been 
constructed. As an extension of studies in the parameter plane, a three 
variable parameter system has been approached by using basic descriptive 
geometry properties and simple illustrative examples have been worked 
out in the parameter space. 

Finally, a real life system with time delay has been analysed and 


redesigned by parameter plane methods. 
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I. BASIC PARAMETER PLANE TECHNIQUES 


A. GENERAL 


Parameter plane techniques are algebraic methods designed to study 
the behavior of linear systems through their characteristic equation, 
when two of their parameters are adjustable. So in fact they are an ex- 
tension of the Root-Locus technique where only one of the parameters of 


, 


the system is adjusted to meet any desired performance requirements. 
P. SILJAK'S APPROACH 


Parameter plane techniques have been at first studied by Mitrovic 
who considered that the parameter variations directly affected two of the 
coefficients of the system's characteristic equation. 

Then Siljak extended Mitrovic's work by including variations on the 
system's parameters "oc" and "R" which appeared linearly in one or more 
of the coefficients of the characteristic equation, which was of the 
form; 

n 
= k 
File) = a,s = 0 αυ 
К=0 ° 


where a, = b ° + су В + а (І-1а) 


k 
with b., c., d. constant and real 
Y i i 
and *,B variable and real. 
Curves obtained in the 4,8 plane are called parameter plane curves 


and they provide information about the roots of the system's character- 


ШОС ециа[їоп іп а similar way to root loci. 





Siljak's work can be summarized in its goals as follows: 

First make an appropriate substitution for s in equation (I-1). This 
substitution must be in terms of C and W» ог О апа Ш, ог {η general in 
terms of quantities one wants to study. 

Second step is to separate the result of equation (I-1) after the 
substitution is made, into real and imaginary parts thereby obtaining 
two equations in @,B, C, w and the constants of the system. 

maty, one solves for ©, B in terms of Z and w. If, for example, 
one fixes C = t and varies U over the range 0 < w < co, one will obtain 
in the o-B plane the necessary relation between œ and B in order to have 
a complex root pair of damping ratio = E (a constant £ curve). ΤΕ, 
on the other hand, one fixes ü = v andgvarics (r over the ranger P s < '1, 
one will obtain a constant EE LUE denoting the necessary @-ß relation 
to provide a complex root pair of natural frequency ü = Wo: Siljak used 
the following substitution in order to study constant and constant ш 


curves: 


nor A-C u C] (1-2) 


where U, and Т. are Chebyshev functions of the first and 


second kind respectively, given by the following recursive relations and 


initial conditions: 


Па 20 +, = 0 


(153) 
бела о ои оед 
where 
T(t) = 1 Ut) =0 
= and of (1-4) 
T, €) = Y U (O =). 


|2 





MENS to be noted that: 


T CO = EDO 
(1-5) 
k+1 
Uu CD = (-1)** u @) 


In the second step after substituting equation (I-2), [considering 


also the results of equations (I-3) to (I-5)], into equation (I-1), Siljak 


obtained: 
n 
Y κ к Е 
συ а 0) U. E) - 0 апа 
Eo 
- (1-6) 
k k 
) (1) аш“ ш @ =0 
к=0 
where a, = ы + c, P + d. (Т-1а) 
So finally: 
ο. =D 
(T-7) 
Вис де spo 
from which: 
тыла ο E Eu 4 оо N 
s вол 2 
where: 
n n 
k k k k 
В, -X (-1) b, W - В, =) (-1) В U 
k=0 k=0 
n n 
к к К =) kayik 
EA C, - ) Cl)ewu. 
k=0 k=0 
n n 
k k k k 
D; =) (-1) d u 1 D, =) (-1) dw UL 
k=0 k=0 


The study for conditions in the a-8 plane giving a real root at 


S = 0, was much easier. In fact, by substituting s = © in equation (I-1), 
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Siljak obtained 
n 
Jead 0 (1-9) 
К=0 


which by using (I-la), led to: 


п п n 

k у E ко 
о) ъс + В се +) де = 0 
k=0 k=0 =0 


k 
which is the equation of a straight line in the cG-ñ plane, for given "Оо", 


Details on this work can be found in reference [1]. 


б. HOLLISTER'S WORK 

In section 4 of reference [2] Hollister obtained parameter plane 
equations for a system with characteristic equation of the form of equa- 
Etom (1-9), with: 


= bÀ + c 


a, = b. в+а 


Ë AB + Λι. (1-10) 


k 
His solution for constant "ш" and "Z" curves was obtained by simul- 


taneously solving two equations of the form: 


PA TG B +D AB + F 0 


ү 1 Е 1 


БА + CB E DAB T E, 


(I-11) 
0 


where Ba, С. D., Е., ӨЗЕН ШІ 21-52 functions of "ш" ана "Са 


well as of the first and second kind Tchebychef functions T, (5 and uU, (D. 


В. PRACTICAL CONSIDERATIONS (AN INDICATION) 

Although the concepts governing analysis and synthesis of feedback 
control systems, once established are quite straightforward, there is 
some difficulty in making an appropriate selection of the range of A and 
B parameter values to be used in the digital computer program. 

Some practical ways for defining a set of approximate ranges of the 
parameter values for use in the first run of a parameter plane computer 


program, are presented in section V. 
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js YSI DO ITL TIME DELAY -- A RESUME 


The purpose of this section is two-fold: 

est: To introduce the reader to the State of the Art of analysis 
for systems with Eransport lag. 

Second: To define the objectives of the work presented in this 
thesis. 

Since stability analysis will be treated separately in a later sec- 


tion, items concerning this subject are not included here. 


A. THE TRANSPORT LAG - AN INTRODUCTION 

l. The existence of a time delay effect associated with any sys- 
tem, linear or non-linear, is inevitable since the working signal (com- 
mand or information), needs a certain period of time to travel through 
and reach the output. So if a signal r(t) is applied in the input of 
one of the most elementary systems, as shown in figure 2-1, the output 
ee) will be in general equal to: 

Е ЕСЕ T) (11-1) 
where T is a time interval associated with the equivalent length of the 
system ''£' and the average speed u of the signal through the device. 
Specifically: 


(II-la) 


= 
N 
s 


There are of course cases when the dimensions of the system and the 
speed of the signal are such to imply a negligible time delay effect as 
for example in the actual lumped A.C. circuitry where the wavelength of 


the signal is large compared with the dimensions of the circuit-system. 


15 





SIMPLE LINEAR 
GAIN SYSTEM 


IN T-DOMAIN 


¢ —_4 


-т---- 


FIGURE £- 


. 16 





But there are also other cases when the time delay effect is import- 
ant and should be taken into consideration. 
In such a case one can define a system in the s-domain,equivalent 


to the one of figure 2-1 by Laplace-transforming equation (II-1): 


C(s) = K (к(өзе 79) - (ке 79 )в(в) 
OUS). . -TS ü 
or πο ; Κο (11-2) 
2. After this presentation one can define the time-delay effect as 


"the physical phenomenon of finite time spent by a signal travelling 
through a finite system." For each combination of signal and system 
there is assiciated a time delay constant 'T'" which denotes the time 


necessary for this travel. 


В. THE EFFECTS OF TRANSPORT LAG ON A SIMPLE LINEAR SYSTEM 

Before the State of the Art in analyzing systems with transport 
lag in the parameter plane is discussed, it is interesting to observe 
the results of this time delay in a simple system by use of conventional 
methods of analysis. 

b. Consider, for example, the system of figure 2-la, with trans- 


fer functions as follows: 


A τν 
78 (5+1) 


a. Open loop: GH 
b. Closed loop: --- = —— 
ΤτΕ 
From the root locus analysis of this system as shown in figure 2-1b, 8 
pair of roots at Z = .707 and Ша 1: derived for K= 5. 

2. In figure 2-2 the same system is shown with the exception that 


a time delay is now present in the forward path and the transfer function 


changes to: 


Ша 
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К =т= (ІІ-3а) 


а. Open loop: GH = SFD e 


=S 
- —— —F ке (11-3Ь) 


b. Closed loop: жен 
5 +5+Ке 


lo 


In trying to apply the root-locus analysis technique one is faced with 
Бисер first difficulty since, 


K 


Επ 
5 (5+1)е 


has an infinite number of roots, as the series expansion of Ε implies. 
93 ee are a number of approximation techniques one can apply to 
overcome this difficulty as mentioned below: 
a. For small T and low frequency, 


"c EST (II-4) 


b. hrom the fact that 





) V 
e - lim (1 m ) 
v- š 
у 
\ 
- A 
and e = lim (5 : ) 
v t2 
V 
one can approximate: 
Set 1 Y я = 
е πο T if v is sufficiently large. In 


1 + 2 


order to have an intuitive feeling of how large v should be for accept- 
able results, note that (l + πο = 2 O producing a deviation of 4.47 
from e = 2.72. So, any v < 4 will produce too crude an approximation 
since (1+ ne Δ, 2,12 = 97, 

ο In the case where both the frequency and T are too large 


. . zu SR р : : 
to neglect in the expansion of e ; it is possible to approximate the 
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time delay term by the ratio of two polynomials in "s" yielding to a cer- 


ШЕНІ order of terms identical with those contained in the series. That 
is, 
N +Ts 
T ) 
D +Ts 
a ) 


-Ts 
е == Е 5) i 


where "a" and 'b" are the order of the numerator and denominator poly- 


nomials. So since 


2 
pos = 1 - Ts А s LI - 
B 2 24 
if third order terms can be neglected, it can be shown by straightforward 
division and a series of comparisons with the series terms that a conven- 


ient approximation is the following: 


Is Даме мо _ 1- 3/⁄4(Ts)+ T /& (e°) 
= l +s ü l + Т/4 (в) 


These are the basic features of the so-called '"Padé Approximation Tech- 
nique." In appendix A the most often used approximations are listed 

(degree of higher order polynomial equal to 12). In this introductory 
study of time delay effects the Padé approximation will be selected as 


the most accurate. From appendix A then, 


Re pas a ` 
“a D (Ts) 
4,4 
1680 - 840Ts + 1801 s* - 20T^s^ + pg 


1680 + 840Ts + 1801282 + 20T?^s? 4 T's? 


and by use of this approximation the closed loop transfer function of 
the system becomes: 
K*N, (198) 


= = — νο. s (11-5) 
Є *s)D, , (Ts) + K*N, д (78) 
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4. The root loci of this system have been plotted for four different 


values of T by use of a digital computer, as the following table indicates: 
TABLE 2-1 


το Figure 2-1b 

T = .1 Figure 2-2a, 2-3a 
T= .2 Figure 2-2b, 2-3b 
ES Figure 2-2c, 2-3c 


T =1.0 Figure 2-2d, 2-3d 


It should be noted that the Padé approximation introduces in general, 
roots in the system which all-together do not disturb seriously the dom- 
inant roots near the origin. Data on the characteristic equation result- 
ШЕМЕСТ each of the above "T's" are shown in Table A2 of Appendix A. 
2. Partial Conclusions 

а. A careful observation of figures 2-1b, 2-2, 2-3, and Table 
2-1 shows some of the effects of the existence of transport lag in the 
linear system under consideration, namely: 

(1) With the insertion of time delay in the system the 
segment of the root ee parallel testhe jw axis doesn't exist, but with 
increasing time delay constant the corresponding segment is bent analog- 
ously toward the jw axis. In other words the system becomes unstable. 

(2) As a consequence the Қ ax of the system which main- 
tains absolute stability is decreased with increasing time delay. 

(3) It is also obvious (Table 2-2), that the relative 


stability is reduced with increasing time delay constant. 
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6. General Conclusions 

From this discussion, the main facts concerning the insertion of 
a time delay in a linear system are apparent: 

а. Апу system carries a time-delay effect with it. Therefore, 
no system can be thought of as absolutely stable. Nevertheless, systems 
with very small time delay deviate so little from the ideal case, that 
they can be thought of as delay- free. 

Б By considering Figures2-3 we see that when approximating 
the time delay term by the Раде technique, poles and zeroes are intro- 
duced at finite W э in the "s" domain. This is true for any other approx- 
imation. 

с. For T increasing, the relative stability of the system 
decreases and the root locus suffers an increasing deformation towards 


БІР ГО direction. 
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C. PARAMETER PLANE METHODS FOR SYSTEMS WITH TIME DELAY (STATE OF THE ART) 


il. The closed loop gain of the system shown in figure 2-4a is 


E ud. 


EC + N(s) 


with n J 
D(s) =) а Е апа N(s) = 55 
k "к 
k=0 k=0 
So the characteristic equation is 
n J n 
Ta x k И NN Каха 
Є ) 98 +) ns = ) ар (5)5 O (11-6) 
k=0 k=0 k=0 


: : Ts à А 
where (s) 15 alfunction of € and Contains the variable parameters A.B 
a, р ; 


of the system. So if only linear combinations of a and b are considered: 


a, (s) = Аз, + ac 9 + Ва + be, e "S + f 4 ge (II-7) 


Pee Int mine to derive constant "C" and "ш, equations for the 


System, one must substitute s in the characteristic equation in terms of 


C and Qs that is by: 
s= -(w + jw / M r* (See figure 2-5). 
n n 
Therefore, it follows (Siljak) that 
* st [s c0 «4 1- uco] 
n k k 


with T, CC) and и СС) being Chebyshev functions of the first and second 


kind defined by 


: -1 
Т, (б) = сов (Ксоѕ Ὁ) and Ὁ (6) = sin(Kcos Ὁ) 


h “ο. C) 


See section I-B for the recursion relations), and note that 


T (C0) - CD m (6) and Uu C0) - CD" (6). 
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E After the substitution for s in equation (II-6) is made, and 
by using the information in the previous paragraph and equation (II-7), 


the following result is obtained after separating reals from imaginaries: 


| AB, (o. C) + вс, (ы ,C) + р, (0,6) + 


з Гав, бо 6) + вс, о 0) + 2,0 ,0| = 0 (11-8) 
where 
n 
ву =) ш | Cr, Gb ae Pe [ C1)ST @)cos9 -( DRA? v esi )] 
k-0 
n ——rF 
Ci =e | ar eh CET, Ecos -ept i-e? u, (©) sind )| 
k=0 
Пп — 
р, =) ® | cor qt s g Í CAT (С) сове НСО еее ЖӘРІЗІ 
ке0 
und ^, 
= уш, “le ιο, o/o ree (a C)/1-¢7cos9+(-1) 7, €) sind | 
k=0 


* i as yde o er coso) sind ) | 


n а 
D, E yw особи од? 7? ( AD cos, ©) sind.) | 
k=0 


sT 


{Note ENe identity € = бш ρον + jsin0) for @ = Cu T апа Ө = о (1-0 r} 


&. The solution for A and B easily follows to be: 


Шо МЕ л лм yc ο; 
o A { Е A ; ши по д 


5. This derivation is the one followed by L. Fisenberg [4]. A 
derivation of parameter plane equations for a system with time delay and 


two variable parameters included in non-linear combinations, will follow. 
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б. System αυ Produc L AB Present in its Characteristic Equation. 


Consider the system in figure 2-4b with closed loop gain 


E s. 
5 E + G(s)H(s) 


For 
Gy (S) H (s) 


G(s) 6, (5) апа Н(5) пъ С) 


the characteristic equation will be: 


Ts 


+ Gy(s)Hy(s) = 


DIS 
E 
o 
o 
[| 
© 


F(s) = Gy (s)Hy (s)e (II-9) 


k 


I 
e 


where: 


Ts 


a, = Alb, + се“ Ὃ * B(d, + ес! °) + ABCh, + те? E + g£ ) 


An expression for 5 can be derived by using 


eo a 


СШ = Ep + JE e ieee 


and calling 
Cw T= 6 and w fre? [у= 
g^ = Ф 5 == 
So 
m Y 51 = e (coso + j sind) 
The expression for бы ав а function of 9» É, and the Chebyshev 


Bemetions are also available. So the characteristic equation can be ex- 


Pressed as a function of Č and Y, as follows: 


N 
μα EE en МЕ Е (Š) - iie U (94 (11-10) 
Den Zen k k 
k=0 
where 
ak = АФ, + сє ©сов®) + B (d, + е, ?со50) + (Еве созд) 
+ AB (h, + ie “cos8) + je "sino (Ac + Ве + ву АВ) 5 antja, 
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Next, by defining PR = т (С) апа Вт ше U G), the characteristic 
equation can be expressed as: 


n 


Fa) = Се)“ ЕЗ s d | 2 
k=0 


which after separating into reals and imaginaries leads to: 


N 
SE [res] = ) C9, GB BER - 9 
k=0 
N Т 
апа Р Im | F(C,w,) | = ) Cu ) (o B - Br) = 0 
k=0 


And finally, by expanding terms inside parenthesis, this set of equations 


can take the form: 


AB, + BC, + D, + ABE, = 0 
AB, + BC, + D, + ABE, = 0 
where: 
N ER Асы. 
B E C E πο e en (C)cosd + U e -e sind ) | 
ще n кк ο. οκ Ben 
k=0 


a—— --- 


N 

Cue )C») à m. (C)+ ee (4, (C)cos® + U, je sind | 
k=0 | 
N 

Di = ) o 08 т (Оз в,67% т. (6) сово + v, «i-e sin) | 
k=0 
N 5 IRE ЕЕ, 

E, = ) Co т (+ “т (О)совд + 0 ОС sine )| 
k=0 


N mee’, Eun 
B, = y a A -ы (61-22 + ce τ. (E) sino - MIES x] 
k=0 
N 


E" ο EO 4 + гє n (sio - v. o Jie 2221 
k=0 
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N ———— — 
D, = ) Cw | - 8 u Ле + ge P(T Csin - v, «i-e cos )| 
k=0 


N 
O) 
E, = )( 994 во (ОД. Cod L € T C)sim UL (E) 2012050 
К=0 
The selection then of (TI-11) is a matter of algebraic manipulation. 


un The derivation in paragraph 6 above is the one followed by 


S.E. Lamberski [5]. 


e. THE OBJECTIVES OF THIS WORK 
1 What has been discussed up to this section was a review for 
people familiar with control theory and hopefully an introduction to the 
state of the art to others. 
ДА The sections to follow contain material associated with the 
fulfillment of objectives as follows: 
а. Theoretical study of general systems with three variable 
parameters and with time delay. This is shown in Section V. 
ον Study of actual engineering problems involving two or 
three variable parameters in systems with time delay. 
er (1) Since for such an engineering study to be complete, 
a stability analysis is needed and hence Section IV is devoted to pre- 
senting the state of the art. 

(2) In Section III parameter plane equations are derived 
for a general system with two variable parameters and time delay. The 
reason that this derivation has been included is that it shows an inter- 
esting way to deal with the time delay difficulty. In all of its other 
aspects it is the same as the methods mentioned in this section, since 


they involve the already familiar algebraic manipulation of Siljak. 
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III. SYSTEMS WITH TWO VARIABLE PARAMETERS AND TIME DELAY 


The purpose of this section is to derive parameter plane equations 
for a general system with time delay and two variable parameters. Тһе 
basic parameter plane method is used, but it is different ‘in the manner 
the time delay effect is treated. This approach will require a shorter 
computer program to implement. 
Since the investigation of this kind of a system started with the 
analysis of a specific actual problem, equations for that kind of a pro- 
blem are first derived. Extension to the general case is then covered 
in two further steps, by: 
first, considering other possible real life system configurations 
and systematically filling in possible terms for a more general 
representation of all of them by one equation, which is true 
under certain conditions, 

and second, removing the imposed conditions and thus getting the 
general parameter plane equation for a system with two variable 
parameters and time delay. 

A. A SYSTEM WITH TIME DELAY IN THE FORWARD PATH, AND THE VARIABLE 

PARAMETERS IN THE FEEDBACK, 

In this section, parameter plane equations for a system with time 
delay in the forward path and Sr variable parameters A,B in the feed- 
back path, will be derived. 

1. The System's Block Diagram 

Assume that the system under consideration has a block diagram 


as in figure 3-1. 
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ШОУ ОГ IIS out that the transfer function C/R is: 


> е8 (в) 


D(s) + a ON + B)N(s) 


хо 


(CE 


Therefore, let As + B = AN, (s) + BN, (s) and assume that the coefficiencts 


n 4 2 
of D(s) are EM щот N, (s) are Legh , and of N- (s) are ο... , 


k=0 


where n is the order of D(s) and £ is the order of N(s), increased by 


one, as N, (s) = 5 • №5). 


D. Characteristic Equation 


It follows from the above discussion that the system's charact- 


- 


eristic equation can be written in the following form: 


n 4, 
y k NV u -Ts k - 
24,5 + ае 5 Ξ0 (III-2) 
к=0 k=0 
where а + Ab, + Be, (111-3) 


D. 





ШОО Л Г ЛОП ОТ a set of Equations That Can Provide Constant ie 
or Constant "P" Curves in the "A B" Plane. 


С ааа 
Starting now from equation (III-2) one can express the complex 


quantity "s" in terms of a and Ç and then attempt to get two equations 


Еле form: 


| 
O 


A*B (Е) + вс Е) t D (o 50) - (111-4) 


| 
O 


A'B, (W 0) + BC, (Ws) + Do (90) = 


One can see that the solution follows Siljak's basic steps. In this case, 
Ейеп, a solution for A and B is possible, for a given O, pair. Easily 
then the constant СЕ and the constant "С!" curves will follow by keeping 
one of the two variables constant, and varying the other over a desired 
range of values. - 

Shown in Appendix B is a way a deriving equations in the form 


пе ones in (III-4). In fact one gets: 
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[у КК + ΠῚ ка ка МЕ PER n. €) 


n o (III-5a) 
4 
(езе a en tp atm) ^ 
2 
4 
JE P" T š M ο ο δα: (с) + (III-5b) 
k=0 


t 


| уа (со + zO = 0 


where Бү: В» c, are constant coefficients of the characteristic equation 


of the system as mentioned before, and T 


kel? kel are given from the 


following recurrence relations: 


iced с т 


(ІІТ-6) 
кк к 
and subject to the following initial conditions: 
= a 2 ; Г | / ú Г ТА 25s 
X = \/1- sin Tu 1-C tC cos Тг (1-6 M 
Em. c 
а = cos; Tu /1 С 
FK. ^ 7η [ το; E 
2 = (А € sin To /1 С | = еск Twy! C j 
ΤΙ = -Ç 
T = 1 
O 
Ti = -ζ (111-7) 


By comparing equations (III-4) and (111-5), 
¿ 
= Та пе ) 
В -е Уьш "X Qu) 
k=0 
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and the rest of the values for LB. ‚С, ΠΝ rl. follow easily. Soi £ 


i=1,2 


В. С, - BAC, # 0, from equation (III-4) follows: 
A= (C,D, - 6251) E (B,C, - В,С) 

(TII-8) 
DNE LO 


4. Derivation of an Equation That Can Provide Constant 'O' Curves. 


The real roots of equation (III-2) may be found by substituting 


s = 0, the required root. It follows then: 


n 2, 
СУ; Jk 
7919 + /е a o O 
k=O 5-9 
4 
where a, = А b. + B C, (III-3) 
Or n / X 
k C -© k \ -© k 
до + де Ab, O + /8 Be, 9 - 
k=0 k=0 k=0 


ce “TO : ; 
and after bringing A,B, and e outside of the summation operator, we get 


finally: 4 2 E 
k 


ИТЕ 
ale a | + В га + jd us = 0 (III-9) 
k=0 k=0 k= μὴ, 


D. Computer Program 


In Appendix D, a computer program able to draw constant O, É 
and w curves or combinations of them is listed. Computations take place 
as follows (Notice that reference is made to equations of Appendix B): 

a. The Case of Constant "Ç" and nw" Curves 


αλ... 6 


аа p 


C and D, are calculated from Equa- 


ү Bo 2 2 


tions (B-15) for the constant quantity of the curve wanted. For each 
; г ues 
value of the running variable, wa 5 5 X1? T. and T.-1 are calculated 


by iteration from the appropriate recursion and initial value formulas 


ДА 





as mentioned before and multiplied by the proper b 


Tu E 


κ. ει апа а, as necessary. 


The quantity e Шол ы саси дгеа to form Ene 


-1 2 
(2) By using equations (B-16) a point P(A,B) is drawn on 
the A,B plane for the constant value curve in discussion, whenever the 
running variable takes one of the values we specify it to span. 
Би The Case of Constant "С" Curves 
As one can see, equation (III-9) represents a straight 


line in the A,B plane for given "0", So in the computer program of 


Appendix D, the quantities, 
А, £ n 


EN kul. k у, К 
e : L122 3 ЕС апа d o 
К=0 k=0 k=0 


are calculated for selected ''σ'', and constant € lines are plotted by 
using appropriate intercepts. 
τ Basic Inputs 

In addition to the types of curves, the selected values 
needed, the span of the running variable and some other information пес- 
essary to o му unit one must input the order of the character- 
ШОО equation (which coincides with the order of D(s) for the most 
usual case when order D(s) > N(s) the maximum degree of the exponent 
in "s" that the time delay term, ex. appears in the characteristic 
equation, the time delay constant "T" and the coefficients of the system 


b and d 


E» Ck k' 
B. OTHER PRACTICAL CONFIGURATIONS - GENERALIZATION 


T. Other Possible Time-Delav and Plant Configurations for Feed- 


Aou i bei dici T YUY E T UY l ыл, Il ан Жаны Аай к с мыны is се алғыс самалы 


Back Compensation (A and B Shown in the Feed-Back Loop) 


In order to get a generalization of the parameter plane equa- 


tions for a system with time delay, one should devise possible configurations 
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of the plant, the compensation unit and the time delay and after solving 
for the parameter plane equations for each case, should then try to for- 
see how these equations would look for the general system with time delay. 
D In Figure 3-2 the first of a set of three configurations 
is shown. But noting that the loop gain of this system is exactly the 
same as for the system discussed in section IIIA, one shouldn't expect 
Eur fferences in the characteristic equations of these two systems. 
Therefore, the characteristic equation is D(s) + N(s)(As + DS cum 
b. The next system to consider is the one shown in Figure 3-3. 


Since for this system: 


N(s) 
D 


5 
а ие) | две 15 τ B| 


D(s 





Le 
R 


4. 


it follows that the characteristic equation for this system will be: 


D(s) + B N(s) + E N (S) = O 
So the basic difference with the original system of section IIIA, exists 
the fact that the term multiplying the variable B is N(s) for the 
Matter case instead of Ne Chae for the former. As a consequence of 


that and going over the procedure of section IIIA, one can see that no 


- Tv 
e Mund term appears in the set of equations (III-5) as far as С, апа с, 


аге concerned. Furthermore, only T U, terms will show up in the solution 


for С, апа E As far as В), В) and D, > D, are concerned, there will be 
no change from the solution of section IITA. 

er The last system to be considered for the feed-back com- 
pensation case with A and B in the feed-back loop, is the one shown in 
Figure 3-4. In this case: 
пралина сә rT 


UON CS Be 5) 


WIA 





BLOCK DIACRAM TWO 
(NO PRODUCTS) 


FIGURE 3-2 





BLOCK DIAGRAM THREE 


EOE ОХЕ 5 ) 


RO), 





MESTRE SS 





BLOCK DIAGRAM FOUR’ 


(NO PRODUCTS) 





FIGURE 3-4 





Therefore, the characteristic equation is: 


i = 0 


"πώεα 
Again by comparison with the characteristic equation of the system dis- 
cussed in Section IIIA, one can see that the difference this time exists 
in the term multiplying the "A" variable parameter. In fact, instead of 
νε”) ie 15 А (в М(з)). So following again the procedure in 

| И -Tw C Е 

Section IIIA, one will get no e n? term appearing in the set of equa- 
prons (Iil-5) as far as B. and B, are concerned this time. Furthermore, 


1 2 


ШОТ U terms will show up in the solution for B, and B,. The solu- 


k 1 2 


tions for Ci» C,» and Di ; р, will remain the same as for the system of 
Section IIIA. 
dr Conclusion 

The solutions for the parameter plane equations for the 
systems described in paragraphs IIIBlb, IIIBlc, remained essentially the 
Same as for the system in section IIIA, with the only difference of minor 
changes on the C. Or В, coefficients (i = 1.2), whenever the time delay 
Peter doesn't show in the characteristic equation to multiply one of 
the two variable parameters A, B. These results are tabulated in Table 


В - 1. 


2. Possible Time Delay and Plant Configurations for the Case of a 
Lead-Lag Compensation on a Fixed Feed-Back System. 


а. Consider first the system as in Figure 3-5. The transfer 


function for this system is: 








(опа (ей -15 
: E + B D(s) = 


= (s+A)*N(s)*e 5 (III-10) 
- 15 aS 
sD(s) + e kes*N(s) + BD(s) + e A(K*N(s)) 


£ 
d. d. ει b. 
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By examining closer the characteristic equation of this system the 
following facts should be noted: 

ΠΠ’ επὲ terms multiplying the A and B variable 
parameters are concerned the situation is similar to the one examined in 
paragraph IIIBlb with the system shown in Figure 3-3. 

(2) For the first time, terms independent of A and B and 
of the form ци x(polynomial in s) are encountered. 

р» In the case where the time delay appears in the feed-back 
path with all the other components of the system being as shown in Figure 
5-5, the loop gain, and so the characteristic equation, remain unchanged. 

En The characteristic equation of the system in Figure 3-5 


can be formulated in a more elegant way as follows: 


n n n n 
k АХ Ts u k V κ. ) οτε. m | 
уа 5 + /е ds + B cs + А )е bs = 0 (III-11) 
k=0 k=0 k=0 k=0 
where d а, E and b. are the constant terms encountered in the follow- 


ing polynomials in "s", respectively: 


sD(s), ksN(s), D(s), and kN(s). 
The fact that the summation operators shown in the equation (III-11) are 
ШЕР арсгагсіпр іп a range from "OQ" to "n" is perfectly legitimate and what 
is meant by that is illustrated as follows and in relation to the system 


of Figure 3-5. 


(1) Consider: k =1 and N(s) É 2s + 3 
ше ДЕ + s. + 65 +7 


where order D(s) > order N(s), as usually happens with most plants. 
(СКЕ О Шоу that: 
2 1 
(a) Site Ξ ms E Ба + 65 + 75 + ος. 


(b) ksN(s)- DN 4 e + ο. 
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A. 


О mee TD ρε sl IM 


1 0 
(d) kN(s) = 2s + 3s, 
where polynomial of max(order) = 4, is sD(s). But for the purpose of 
neater summation operations in a digital computer one could form the 


characteristic equation as implied by equation (III-11), as follows: 


n=4 n=4 Е“ n=4 
= N = “Ts, е k \'-Ts, К 
) = 2 
2% ES + де je + B ας + А де Ь 5 0 (III-lla) 
k=0 k=0 k=0 Е- 


мһеге E di С 2 b. for the system are tabulated in Table 3-2, and n=4 
is implied by the fact that the polynomial of maximum order for this case i4 


e Order 4. 


TABLE 3-2 
Polynomial k 

4 3 2 1 0 
5) (в) d, = 4 d, = 3 1, = 0б di = 7 d. = 0 

киы: 1 = ! = 1 = 1 = 
k sN(s) d, = 0 dz 0 4, 299 3 d. 0 
D(s) SE 0 C= 4 Cy = J с, = 6 m 7 
kN(s) b, = 0 b. = 0 b, - 0 b, = 2 b. = 3 


3 Generalization 
ас Having in mind the form of equation (III-11) as well as 
the results tabulated in Table 3-1, one can deduce that a more general 
form can be given to the characteristic equations of all systems discussed 


up to now in section III and this is as follows: 


n n 
Т -ба2*Т*5 "T: je k a op Tes, К 
ж Ὁ ds + В )е cs + А je 5,5 — 
E κ-ο k=0 Aa 


where d с. апар 


κ. αν» εκ are defined as in the previous illustrative ex- 


k 


ample of this section (p. ), and EE δ, δ. are variables that саке оп 


values of 0 or 1 depending on the system under examination. 


5% 








D. Constant “Ὁ Curves 
Formen ЫШ Опо gets by direct substitution from (III-12) 


the following equation representing as usual, a family of straight lines: 


n n П 
к-0 к=0 ке 


с. Constant αν or Constant "ш" Curves 


It has already been defined that: 


0 
ο... - l or 
i i(i=d,,c,b) 1 
ШЕ now, о; Бает пеа аз . 
6! | Қ 
ου. 
Ша то Томиз that: 
Ome On = 1. (III-14) 
1 J 


So by using Š, and ò; one can selectively maintain or discard the differ- 
ent terms that must or must not show up in the parameter plane solution 
of a particular system if this system has a characteristic equation of 
Шіс storm as іп (III-12). 

(1) Example. By examining the solutions for В, апа В, 


in table 3-1, it follows that for the systems discussed up to пом, В, 


paces one of the forms: 


4 
(а) В = уь ш К 
ШЕ Сат E 
k=0 4 
or (b) В, = за реа 
k=0 


This statement can be reduced into one single equation by using (III-14), 
ead tne fact that В, ттге to b, S and therefore to the equation 


ор =] (observe δ. to be related with "bi coefficients in (TII-12). 
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n n 
-0, * T*Un*C N. k k 
= b DES x E 
(ο) В, е Š, Ab Σι δω. T. (ІІІ-15а) 
k=0 I0 
(d) Similarly: е 
-Ob*Tewy °C N k 
DONNE ос отар 
n τ. (III-15b) 
x k 
Li 
+ M b (ET, + T, 4? 


k=0 


Solutions for Cy > с, апа D, > D, follow similar patterns. 


ССО interest to note that if d. = 1 (57-0), 
the presence of X's and X, 1/5 is implied in the relevant expression 
{ IE m = 
for В, B,» С.) 022 ΟΥ ру, р. een Ne Contrary, 9: 1 (б. 0), the 


1 1 . ° а 
presence ofT s and Ті s is implied. 


(2) Derivation of Equations 
For one to actually solve for the parameter plane 


equations, one must make use of the fact (as in section IIIA), that: 
в + w (cosd + jsin8) = w ak 

n n 

and apply it in equation (III-12) to get: 


n n 


y k "d аи "бо" Те (созб + 151п0)00) κ, 1κθ 
даа, (cos kð + j sin КО) + дау (е Io e + 
k=0 k=0 
n 
(ве 91со89 + jsi) ). " k ike 
k n 
k=0 
n 
τα = D | > я 
(де Op T(cosd + jsin® on), " k ike 
k n 
к-0 


which, with appropriate manipulation inside the summation operators in- 


+ 


= 0 


volving dy c, and b >» yields the result shown in equation (111-16). 


(See next page.) 
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| 
I 
| 





п 
К 
i sin k 
ya, (cos k9 + j sin КӨ) + 
k=0 


'Dd2* TUE ΣΝ ἴδω" cos(k9 + wT) + O )| + 


k=0 
Ἢ ' ] 
j ^Ш sin(k9 + аме + ба вт (КО) |) T 
DE D. ES t vk fo, ° cos (k9 + ου + O: ‘соз (к9) | А 
E 
j B t sin(kð +w TE) +6! sin (ko) |} 
n 
OH? Tew f | : т 
Ае b E E а / Py, k δ. ° cos (k9 + Το + M . Ja + 
k=0 


j |5, ° sin(k9 + port + δι sina) ]! | (111-16) 


From this stage by using appropriate substitutions for cosines and sines 


with arguments (k9) ог (К + хит"), in terms ofX,'s,Y,'s, ET and 


U 's and their recurrence relations, exactly as defined derived and 


k 


used in section IIIA, one can form the following equations after separ- 


ating reals from imaginaries: 
n n 
k δη Tem. ЛОВА k 
W Z . 
24, ШЫН ΤΟΠ. 
k=0 | к=0 


E 


n 


© Ten τ \ k 
C n X ' 
Be M (δ. -— от.) + 


k=0 


peb Tn E Sp o K бх +6т) -0 (111174) 
k=0 


(see next page) 
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апа: 


M k да?" Теш “Е 57 Т k ! 
Ши щете ΤΌ 


or tau aa 


k ' 
Be UJ (0 Y. + Ó U) + 


s 


Корт С до Т co би ьш) = 0 (III-17b) 


Equation (III-17b) can be expressed in terms of k Š and T. S instead of 


U's and „55 in exactly the same manner as in section IIIA by using 


equations (B-9c) and (B-12c) of Appendix B. So the second equation of 


the parameter plane system becomes: 


n 


n 
E \ k Одо" Теш Да tray 1 
ee + E, ste τ γάμο δη, СЖ ua IE a А 
k=0 k=0 


n 
BET Te Ts gx + x.) + OLEH + Ty) + 
k=0 
дед . 
Бет" С Е. " p 4 τ. id "ο" er 4 T. М 200) СІН с) 
k=0 


For easier reference, recurrence equation as well as special symbol re- 


mations are listed below: 


(а) Т, 


ee 


(с) Initial condition relations usually necessary 


-2CT, = 1.1 


m2 digital computer calculations: 
(i) ΤΊ z Ee To = us Ti = ae 


Gi) x, =-| 1-6” іп(т 1-6 сов (то 1-0) | 


X = сов(Ть/1-С°) 


58 





ОЕ Jit sin cm ich - eos (au 1-62) 
(iii) DE Шс! 
Ela 
c { : 


ER jr (d,,0,b) 


(3) Parameter Plane Curves. Constant MD and/or con- 
stant "C"' curves can be drawn in the parameter plane in exactly the same 
manner as discussed in section I-B (Siljak's approach). 

(ЕЕ Conclusions. 

(a) By comparison of equations (III-17) and the sol- 
utions for Β1261951 9500292 tabulated in Table 3-1, one can observe that 
the derived equations (III-17) completely cover each case and furthermore 
they cover the case of the system shown in Figure 3-5. 

(yma inc equations ή). are valid for the case 
that the terms encountered in the characteristic equations examined, and 
involving the variable parameters A and B are either of the form: 


UNS S 


)* (Polynomial, (s))+ B+ (Polynomial, (s)) 
Sr Oi the form: 
A* (Polynomial, (s)) = (Be “°)+ (Polynomial, (s)) 


This was because of the specific features of the systems examined. But 


since the terms not involving A and B in equation (III-11) are of the 


form: 


e^ ^5. (Polynomial, (s)) + Polynomial, (s) ; 


the possibility exists for a more general case with A,B as well as AB 
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Mepearing in the characteristic equation: 


F(s) + гБ + 4 B) + ве (ве 19 | ἡ ο(ϱ) + se Pl, 
' - 15 ! 
АВ D(s) + D' (s)e | - 0 (111-18) 
The parameter plane analysis of such a general system follows next. 


Б. NON-LINEAR VARIABLE PARAMETER PRODUCTS INCLUDED IN THE CHARACTER- 
ISTIC EQUATION 
E In figure 3-6 a two variable parameter system is shown, which 


will result in a characteristic equation, containing the AB product of 


these two parameters, as follows: 


F(s) = | D(s) ч sN(s)e 19 + A М(в)е 7] + B s^u(s)e 1°] + AB sN(s)e 7* | 
= 0 
or, 
n+1 n n+2 n+l 
a k k -Ts k. k "Ts k -Ts > k -Ts 
F = с 1 | “с 1 1 
(s) A Cf s + fs e ) + A (by e ) + В NASE e ) + AB / (418 е ) 
k=0 k=0 k=0 k=0 
= 0 { 1111) 
А. However, in order to gain more generality so as to apply the 


pesults into a broader system classification, it is wise to solve for the 


parameter plane equations of the following characteristic equation: 


n- 2 
ет еке са N -Ts. k 
Е 1 ) 1 

)« B E Ее )s + ЈАС + be ys + 

=0 k=0 

m P 

hx Е m я 

/B(e, + cle I + JAB(d, + die Is 0 (177-209) 
k=0 k=0 


and by arbitrarily choosing п = max(n, £, m, p) as well as by filling 


non-existent coefficients of lower degree polynomials by zeros: 


COT 





BLOCK DIAGRAM "SIX - 
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КО), NO | Em (65). 
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-ν Bs 
БЕ 


y 





κα NG) | ‹ } A 





FIGURE 3-6 





п n 
N TS К ADS τ 
Í jí уф p 
ri k + ке )5 + А )( κ t bre )5 + 
k=0 k=0 
n n 
E (7 Ts, k D D ac Ко 
B 2. + се За + АВ) (А + Фе s = O 
k=0 k=0 
j. The solution for constant o and € curves is again based on the 
equations: 


19 
= 8 1 1 = } J 
5 w (cos + j sin 9) ше 


and x= cos(kO + Y TÉ), Mz sin(k9 το TC) 
T= сов КО 5 U= біп КӨ 


(escorted by their recursion relations), being substituted in equation 


(III-20), and thereafter on a separation of reals and imaginaries. 


4. Actually, the set of the two equations is: 
a = = 
АВ, ar BO, + ABD, E, в 
АВ, + BC, + ABD, = -F3 
where: 
n n 
V k τοι k 
= УЬ п ' 
E е μα l 
k=0 k=0 
n n 
x k To P Y k 
= \ QJ nd "ы 
ου е τῳ A 
к-0 к-0 
п п 
ү k ΠΣ. k 
D. = W n τω 
p^, T, em aem x, 
k=0 k=0 
n n 
E k τωρζ y gt К 
35 7 ү ШЕ DAP "E. x 
k=0 k=0 
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n 
k 
в, = Бю, €T, + ioe? η 
k=0 


n 


æ os 


2 
k=0 


n 


— 


mE k т С 
= ; n 
р, guns (CT Em = + е 


k=0 


n 


k 
= W qu 
C C. τ (CT, T _ 


1 


Tu а ү. 
)+ e n ЖЕ” 


К - 
σαι + Х 1) 
к=0 

|@1 
Ὅν... 
Lkn k K-i 
k=0 


n 


— ` k ша ШТ Bray k 
E- EA n CT t T1) +e LI CX, * XV) 


k=0 


k=0 


5. According to Hollister, [2], then: 


= A = 
р -ίδερ + Ac) 8 d py t Ac ep en 


and after B is defined: 





В, E D,B 
where: 
C 
1 
боп D 
1 


2 
2b 
2 
бер 
F, + 
Ba + Ρο} 
Ε F 
1 2 
5 AFD - 3 etc. 
ү ο 
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IV. STABILITY ANALYSIS IN THE PARAMETER PLANE 


In this section the stability concepts in the КОП ОТ plane are 
Ши ТЕ very qualitatively, at first for a delay-free system. These con- 
cepts are based on ideas similar to the ones used in conformal mapping. 

For the case of systems with time delay, Karmarkar's method, based 
on Cauchy's principle of Argument or Michailov's test and its extension 
by Pontryagin, is followed and used to examine the stability region of 


the system shown in figure 6-1. 


А. SYSTEMS WITH NO TRANSPORT LAG 
ІК Given a certain polynomial CE(s) in "s", and of nth degree. Іп 
general out of the n roots of this polynomial, r. will lic in the shaded 


region of figure 4-la and En in the unshaded one, but 


n = ra + E (IV) 
So the entire 's" space is divided into two sub-spaces, by the s = ji 
line and the © = O point which in a sense are their complex and real root 
boundaries. 
2 Now considering the closed contour c, Sconsisting of сі апа 


Бе the jw axis) as shown in figure 4-1аф, we have to realize that it 
maps in the AB space by appropriate mapping of its boundary. But since 
the CE(s) polynomial serves as the mapping function, it follows that for 
different characteristic equations the 0 = 0 and s = jV boundaries will 
map differently. It should also be noted that these boundaries will 
generally decompose the whole AB space into a finite number Ó of regions 


К., ο ο οὗ which will contain í roots in it. So: 


б 


— 
% 


Zi =n ^ (IV-2) 
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FIGURE 4-1b 


65 





>. The whole ЕЯ of the stability analysis is to define the 
number of roots existing in each of the regions R.. o ocne words; "to 
determine the region index") 

а. Contour Shading (The Shading Rule). 

Consider the mapping process of figure 4-lc, where the complex 
root boundary is mapped in the AB plane. So if the mapping device lies 
ME point M, it follows that CE(s) has a coste = ju... In other words, 
maps into "M" via the mapping function CE(s) = O. 

Considering now two points in the M neighborhood, which lie in 
either side of the É = O boundary, one realizes that while moving from 
point M, to М), at least one complex root in the s plane, moves from 


[ 
5, to s, or vice versa. It is therefore necessary to orient the boundary 
Bons, according to the following rules: 
a. m the s plancia right handed rectangular system (t,n) 


is established, where t is the positive tangent direction and n is the 


EE со it in the appropriate direction. The boundary, then, is shaded 


always in the dircction of positive n. (See figure 4-1c, upper). 

Since the complex root boundary in the s-plane is s = jw, 
positive t is determined by direction of increase of the variable working 
along it. 

b. In the AB plane, although the (t,n) system is always right- 
handed in the s-plane, this is not the case for the AB plane. So depend- 
ing on the particular system under examination Ehe positive n direction 
can be either side of t. The positive t direction is determined as before. 
БО гог: 

(1) (t,n) Right-Handed; The positive "n" direction is to 
the left of positive "t", and the boundary is shaded to the left of pos- 


K ve МЕН. 
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ROOT BOUNDARY SHADING 
IN THE “Sano A,B” PLANES... 
— Yi ζω 





р, (-0 ες 
χι. PorT?0) 
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FIGURE 4-1c 
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(2) AI e cs tHemded: Since S "n" will be direct- 
EN Che right of positive "t", the boundary is shaded to the right of 
Positive "t", 
Ec. The Significance of Shading. 
mec haded portions of the regions in both the s and AB 
domains correspond, in the specific case examined, to the side of the 
domain where all the complex roots exist. That is to the stable region 
since the boundary in "s" was the absolute stability boundary. 
d. Rule to Define the (t,n) Co-ordinate System's Orientation 
(1) Consider the polynomial 
Е(5) - СЕ(5) (Ту-За) 
(2) Substitute s - O 4 jU, to get 
Е(с,ш) = R(0,w) +31 (0,0) (IV-3b) 


(3) After forming the Jacobian 


- E ΠΤΙ, с. 
k κ “АЗАЛӘВ/ 7 Әв/әл/ οσο 


decide as follows: 
O (ға) is right-handed, therefore 
shade to the left of positive "t", 
ОБОО ОССО) їз left-handed, therefore 
shade to the right of positive "t", 
e. Similar rules apply for the real root boundary, (the 
O = 0 curve for this case). However, the following rule applies success- 
fully: 
(Е) Since at the point of intersection of s = 9 with the 
S = jW curve W changes sign, so does the Jacobian J. Therefore, in the 
AB plane the G = O curve must change direction of shading at the point of 


intersection with the ¢ = 0 curve. 
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(2) In the stability region neighborhood of the image o” 
gn the s-domain, the shading of the 9 = 0 and € = O curves must coincide, 
Since this is the case for the s-domain shading also. 

5. The stability region in the AB plane, for a delay-free system, 


EEUU been established. For more details one could refer to reference 


mE. 


р. ENSTEMS WITH TIME DELAY 

L, Ши Ерис subsection, Karmarkar's approach, [6], to define if the 
domain of maximum number of roots is also the stability domain, will be 
followed since it provides in a single plot the stability domain as well 
as the test. | 

D. Karmarkar's stability theorem states the following: 


απ τας СЕКС) (С ре а transcendental polynomial: 
m,n 
В 2 р „05 - 
ο 5) Ne e (IV-4) 
P Ir 


where: 


a” are rcal coefficients 


an - Ὁ (πε principal torm) 


p,n are intergers, and 


N τος О 1 ү 05 р 
ENS e = E слага а, ау... аң, А (IV-5) 
τυ i Ain} | ° 
1 Y "ne 
“m0 m2 mn N 4 


For CE(s) to be stable it is necessary and sufficient to find an interger 
K so that, starting with the line B = bo: the total number "A" of alter- 


nating intersection of the € = 0 curve with the A = a, and B = b. lines 
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15: 

A = 4kn + m (IV-6a) 
when ü) ranges over the interval 

ОИ ОО ας 6/2.) for 0 S € < п/2, 
Bere, 050ш52г(к - 6/Л) for 0 < ó <l, (IV-6b) 


In figure (4-2), the results of this analysis for the system 


2 
СЕ (5) = DE + s)e” +as+b=0 
is shown. For this particular system one has m = 3 and n = 1. So since 


Ши was varied in the range ὮΙ = (small + 14) and therefore 
7; Ξε (where w in rad/sec) 
k=2 max 
for k=2, the intersection predicted by the theorem should be 


P= | (кп +) |, = 11, 


Bl 
Ше 


and this is what is measured, so the system is stable in the domain s 
bounded by the © = O line, which is the B=0 line as easily one sees from 
CE(s), and the first encirclement of the (=0 curve since the Jacobian 
Ш 53), is less than zero for this case. 

В. The actual derivation of Karmarkar's theorem is based as already 
mentioned in the following theorems: 

I a lovis Criterion, as a direct application of the 
Cauchy argument principle [71: 

b. Pontryagin Theorem on the zeros of elementary transcenden- 
tal functions,l8], and the conditions imposed on A and B when the mapping 
operating point scribes the ) = O curve. (That is equivalent to having 
roots of the form s = jw in the s-plane). For details on his derivation 


see reference [6]. 
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CONSTANT $:0 AND 6-0 CURVES 


— MALA crm r 


FOR A SYSTEM WITH J DELAY 





FIGURE 4-2 





С. APPLICATION ON AN ACTUAL SYSTEM 
j|. Consider the stability of the system shown in figure 6-1 , for 
different values of the time delay constant T = 0.0, 0.1, 0.2, 032 

It is obvious that for T = 0.0 one can apply methods discussed 
in section IV A, and for the rest of the time delay constants the methods 
of section IV B are necessary. 

However, a far more practical method of determining the stability 
of the system is to arbitrarily pick two points in either of the smallest 
two areas defined by the closed constant € = ONE curve and Ше с = 0 
straight line, and simulate the system in "t"! domain. If for one point 


the system is stable, the domain it represents is the stable region. 


2. Тһе characteristic equation of the system is: 


EC - B 4 98.4696s? 4 4941.1680s^ 4 8978.082s^ 4 5086.9961s^ 4 


7 600115 - 147.59 | onus ES + ES + .488s) + 


ПЕ ο ος + 488) | = (5) + Ba E h(s) | 0 


a. In order to obtain constant "o" curves, let s =O., 


= 


ο ο ος Е g(c) + B πίσ)| - 0 


ad 


and for the O = QO curve one gets: 


СО ОТ ЕШ О = го), 
Ок ООЛО p x (488) - 147.59 


.. for finite time delay, the constant "O = 0' curve is: 


B = 302.87 (VES) 





* 
: Notice the equivalence of the [(=0 and 2 = O curve. In general 


Belt curve is the loci for roots "s - + jw." 
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Б. After multiplying equation (IV-7), by a one gets: 


IS Ces + eS | 


a ^ 


A g(s) + в h(s) | = 0 (1V-7a) 


ЖИН СЕ Т 15 пос ап integer, call T = t,/tg, where "tp" and "tg" are 


integers,. and substitute into (IV-7a): 


! ! 
ens f(s!') + PS É g(s') +B his‘) | ο (ТУ-9а) 
where 2 
r _ 8 = 
s = т (IV-9b) 


As a consequence of equation (IV-5), equation (IV-9a) can 


also be written as: 


(в! ΧΕ E Е ХА + .4888 | nu |- 147. 59 |} n 

(в! Че 5. E 488A + аз | + еба 3162.6011 |} | 

в") 40 [адад + в ο. 1 | 5086. 9961 |} + 

ο ο + 0 x p| 2 P _8978..082 | + 

(ο!) οὔ" | οἱ + οὐ | 4941.1680 |} + 

(te^? | οἱ + e^? [ 98.4696 ]) F 

D a + etn? | ТЯ o]; - (IV-9c) 


According to the definition of section IV B, the principal 


m ns | 
Em (sa e `), is: 


mn 

ο ο οσο (ТУ-10) 
Ете п - 6, п- с anda = 1.0 

SOC is) to be stable 16€ is needed to find an integer 
k such that for w' = (»/tq) in the range 0 = w' = опк+ > = п(2к+ 2), 
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where: 


ОО осе, 


5 


Ku I] number of intersections, A, of the € = O curve with the A,B 


axes moved to a point Р (Ата > В - b ) is: 


À = 4kn + m E (^t )k + 6 


(eee ee ee ee a 


In Table 4-1. the columns under the heading "calculation", 
epee filled in" according to these two formulas giving "X and w" for 
= - N35. which correspond to time delay constant values of T= (.1, 
2). The rest of the columns are filled in relation to a computer 
solution of the "£=0" curve vs frequency and this will be discussed later 
in this section. 

es Another fact of practical importance to note is that 

karmarkar's theorem starts counting crossections of the C=O curve with 
the axes through a testing point in the A,B space, from w=0. But any 
ПИО БАІ machine must have an initial value of #0. 

So it is necessary to calculate the origin of the £=0 
Ее for W tending to zero. This limiting point can be calculated 
rather easily since for the case of T=0 the time delay effect is not 
Present. But also for w > 0 on the (=0 curve, the time delay effect is 


also unimportant since: 


lim ae = lim Eum (IV-11) 


τοῦ WO 
So the origin of the constant €=0 curves is the same for a system with a- 
variable time-delay constant, and independent of T. 
е. For the particular system of equation (IV-7), one takes 
Ше limit by operating on CE(s) as follows: 
limit of ((=0, curve) = lim CE(jw) 
WO w-0 
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his operation will strike out all the "5" dependent terms of CE(s), 
and. therefore 


lim CE(jw) = -147.59 + .488B 
WO 


and by equating that to zero the solution is again 
В = 302.87 (IV-12) 
One should anticipate this continuity of real and complex root boundary 
since for WO the time delay effect is not present and therefore the dis- 
cussion of section IVA, is valid. 
Í. It remains now to determine the value of the Jacobian in 
Pia tion (1V-3c). This Jacobian is in general a function of O, ü, A, B 


ESXror a particular о 


lim J(OWW,A,B) - lim 


бза сета 

К R M UE | R | n 
Lim 97 . lim a — Lim — . lim E 
η б-а 9 gma 9 Oa pr 


and because of the nature of the derivative operation, and the well be- 
haved characteristic equation, the order of the limit and derivative 


operators can be interchanged. So for O-a = 0: 


OR(0,u) CUA), IR(0,w) aT ο) 
ЖЕТЕР 5:5 57 2 


Now referring to equation (IV-3b), R and I will be computed 
for T=0, and the same orientation of the (t,n) coordinate system will be 
Essumed for T = O since for WwW = w, = O and a differential change in W, 


0 
o cause W = Wo + ау = 0 + dw), the (t,n) orientation shouldn't change. 
But once established the orientation of the €=0 curve cannot change. 50 
the orientation for a system with variable time delay is independent of 


Т. 


Now taking into account that some terms of R and I will 
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disappear in the derivative operations, define: 


ОО - - - WB + .488B 
mod 


It 


(0 ,ш) ΠΠ + .488YA + .4440B 


mod 
to be regular R,1 polynomials, with the A,B independent terms deleted, 
amd the expression for J is: 


J(0,w) = ES + ‚290? + 488,2 - .238 


Considering cases with w ^ O: 


lim J(0,w) = у = 
? (roo 

and O O 73909 0 
ur0 


So at the 'w'” extremes J < О and due to the nature of the 
problem the same must occur for intermediate values. Therefore, (t,n), 
is Left-Handed and for w increasing the shading goes to the right of the 
C=0 curve. 
3. Digital Computer Stability Analysis 
a. With all this information available the "o = O" and "£ = 0" 
curves are solved with the aid of a digital computer and the program of 
Appendix D. 
b. timeswasedone for = (0.0, 0.1, 0.2, 0.3), and the re- 
sults are displayed as follows: 
(1) Maximum root regions: In figures from 4-3 to 4-6. 
(2) Stability analysis: For a роп рп Ее А,Б plane 
which lies in the general area of the system of section VI is going to 
to operate, (that is, P = (A = 5722, B = 10568)), the É = O curve was 
drawn resulting to crossections with the A, = 5722 and B, = 10568 lines 


0 0 


as shown in figures 4-7,4-8,4-9. 
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С; By using this information the rest of the table 4-1 was 


filled in and the partial conclusion was: 


(1) For T = 0.1 and 0.2, the system was found to be stable. 


(2) For T = 0.3, the "w range" of the simulation obtained 
showed that for k=0, k=l the requirement for stability was not obtained. 

As General Conclusions 

The first thing to notice is that for T = 0.0 one gets the big- 
gest maximum root region. This AB region is with no other reservation, 
the stable region in the AB plane. 

As the time delay constant is increased the maximum root region 
is decreased. But still, this maximum root region is not the stable 
Eon. According to Karmarkar's criterion the system with T = 0.3 is 
not defined stable for k=l and so one should try k> |. These results 
will be checked out and evaluated in section VI, where this practical de- 
sign problem is encountered in more detail. 

b Notice also that for the one variable parameter, "k", system of 
figure 2-2 the conclusion vas that for increasing time delay constant the 
maximum k for stability was decreasing. So in a sense the stable "para- 
meter line" was decreasing. 

In the two variable parameter system one has also a decreasing 
area of stability in the parameter plane with increasing T. This should 
be expected since the physical situation remains the same. 

As a result of that we should expect the volume or the hyper- 
volume of stability, if the system can be stable at all for a given time 
delay constant for a third or higher order variable parameter system, 


to decrease with "T" increasing. 


ЙТ 





Delay 
Constants 


072 


0.3 


10 


10 


10 


TABLE 4-1 











k Calculation Simulation 
selected À W À W 
max 
3 18 62.8 (45) 18 203.5 
4 
= гад /зес 
20492 
rad/sec 
[| 14 η. 14 yiri. 
7 
= rad/sec 
78552 
rad/sec 
1 18 62.8(1+2)119 77.2 
4 
= rad/sec 
Woe, OD 
rad/sec 


Remarks 


Stable 


Stable 


Not defined, 
Should 
check for 
k > 1 
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V. SYSTEMS WITH THREE VARIABLE PARAMETERS 


The purpose of this section is to provide a general way of analysis 
of systems with three variable parameters to meet specified conditions. 

In the analysis to follow, a certain knowledge of the basis of de- 
scriptive geometry is needed, so for the reader who is not familiar with 
the subject, reference L9] is recommended. However, the elementary con- 
cepts needed later are presented as well. 

The attack to the solution of the most general case of a system with 
time delay is done indirectly. In fact, a system with characteristic 
equation linear in A,B, is first analysed so that the concepts of the 
solution are being built first for this simpler case. The transition 


then to the basic objective is not difficult. 


А. ΙΝΕ OF DESCRIPTIVE GEOMETRY 


E Consider a three dimensional orthogonal coordinate system, with 
ахез А, В, Г. Shown in figure 5-1 is the set of points Р,, Py» P4» Р,. 
2. The purpose of the descriptive geometry is to provide means of 


accurate representation d measurements of a three dimensional space in 
a two dimensional one. So if for this purpose the plane OBI is selected 
as the reference where the measurements will take place, the projections 
of any point of the actual system on the actual OAU plane, can be shown 
on the reference without distortion by just a 90° clockwise rotation 
around the Ol axis. Figure 5-2 then will result in which B and -A as 
well as A and -B axes coincide. The projections of any point on the OBI 


plane are shown by "+" and on the OAl’ plane are shown by "p". 
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2. Points on the Reference Plane 
F poeciioriroliiiSures 5-1 and 5-2 one can note the following: 
ar A point P in the actual coordinate system, can be repre- 
sented by two points P' and Р" on the reference plane. So we have a one 
into two forward mapping ( and a two to one reverse). 
Br ШОЕ ео роіпев (Р.Р!) that form the image of P, lie on 
Ane normal to the Ol’ axis, due to the nature of the image generation 
process. 
Gu If a three dimensional point lies: 
(1) on the OB' plane, its P" must lie on the Ql axis in the 
descriptive geometry representation (e.g. Р,). 
(2) on the OA plane, its P' must lie on the Ol axis also 
(e.g. Ро). 
(3) on the A=3 or A=-B plane, its P' and P' must have the 
same distance from the OF axis (e.g. Ра). 
(>on the Ol axis, its P. and P" coincide (e.g. Р,). 
G. Straight Lines on the Reference Plane 
If now points Р, Р, аге joined together to form the straight 
line they define in the actual space, it is reasonable to anticipate that 
the images of this line in the reference plane must be the two straight 
lines that are defined by (Р,' Pz") and (P," P4"). Furthermore for each 


pornt Px ) on the three dimensional space there is one and only one 


(а,р,с 
Bet of points Px' and Px" in the reference plane such that: 

a. Px' lies on the P, Pz line 

be Px'' lies on the PUP line 

ος If g is the intersection of OG and P' Px' on the reference 


plane; 


lg Px'| ΞΡ. lg psu =a, |0“ 





E 
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о. Planes on ο ae once Р1апе 

A plane is represented by two of its lines or three of its points 
Ш (Пе 3-D space. So it needs four straight lines or six points in order 
to be defined in the reference plane. Notice the ee when the 
actual plane is: 

a. Normal to the Ol! axis. In this case all the four images 
of any two lines selected to represent it, will coincide and will be nor- 
Mal to the Ol axis. 

b. Represented by its intercept lines on the OBI and OAT 
planes. In this case two straight line images will suffice and the 
other two will coincide with the Οἱ axis. 

6. Lines Other Than Straight, on the Reference Plane 

Much of the power of the descriptive geometry representation is 
lost when one deviates from straight lines or planes. For e 
intersection of a straight line and a plane can be determined by a 
standard procedure, but of a curved line and a plane the intersection 
can be determined by trial and error only. 

However, any continuous line in actual space will give contin- 
uous images in the reference plane and vice versa. 

So if one knows that the images (P,', P1"), (P, > P"), (Ра) P4") 
Pas and 


Р, which formed a continuous line in space, then one can anticipate that 


and (P, ыр,” as shown іп figure 5-3b projected by points Ру, Р, 
the projections of the actual line segment would more or less look as 
shown in figure 5-3b and the approximation of the actual curve shown in 
figure 5-3a would be more accurate if more points at closer distances were 


available. 
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В. GENERAL THREE VARIABLE PARAMETER SYSTEM WITH TIME DELAY--A CONCEPT 


ү, A simple extension of the case of a system with two parameters, 
entails that for the case of three variable parameters A (Alpha), B (Beta), 


I (Gamma), the general characteristic equation would look like: 
| - Ts 1 H - Ts i I - Ts | 
4 C, (s) + е С. (в) | + D, (s) + е р, (в) | + E (s) + e E(s)| + 
А Р (5) + me (s) | + A H (s) + EE (s) | + uil x (s) + брі (s) | + 
1 2 iv 2 Ex 2 


ABI] K, (s) E e T^y, (s) | Т |, (ο) + BL, (s) | = 0 (V-1) 
Да The solution for the ČC and Ша curves for this system must follow 
the following steps: 
a. First derive two sets of equations (separating once again 


reals and imaginaries) in terms of € and Woy eg. 

VU Š ») + Ге J q 1 n Ty Io 4 
Ax, (5,9) * B5, Q,5 ) + Ге (Сл) + АВФ (С) + АЛ Кл.) + Bie Qr) + 
АВГр ] (С 0.) = | (С ‚ш Ὁ) > and (V-2a) 

"n : А \ η 27 rey 4) y 4) 
Aa (62) + в6 Фр.) + Ге (Со) + два, Фл) + Мб) + Ве Фо.) + 
АВ (Сю) = -5,(6,0 ) (V-2b) 


b. Follow che same steps as in the two parameter case by: 
(1) Giving actual values to W and Ọ to form the certain 
Seonstamt performance characteristic curve," therefore making αν]... 
апа Yo о... constant for each "бш" рали 
(2) Solving the two equations of the form of (V-2) but 
with constant coefficients, for the three unknowns A, B, Г, Бу applying 
some restriction on the solution, e.g. by watching a certain plane cut of 
ШО А, В, Г space. 
сі Put the partial plots together to reconstruct if possible 


the whole picture in the A, B, T space. 
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B. This idea of taking cuts on the three parameter space is not new 
and in fact solutions by quantizing one of the three parameters (which is 
equivalent to subsequent cuts normal to one of the three parameter axes), 
have been successfully presented. The only difference is that this 
character of analysis has a trial and error characteristic which of course 
is reduced if the physical limitations of one of the variables are more 
or less narrow. 

4. However, by taking enough cuts with planes passing through one 
of the coordinate axes of the parameter space, one can solve an equi- 
valent number of two dimensional problems which can be displayed on the 
same reference plane by use of simple descriptive geometry ideas. 

Э. These ideas can be very easily grasped for the case of a sys- 


tem simpler than the general. 


С. A SYSTEM WITI LINEAR CHARACTERISTIC EQUATION IN A, B, AND I 
- Consider, now, equations (V-2) for the special case that AB, 


fees, АВГ ао пос exist. In this case: 


де, ув у + Bd, (ув ) + Te, Cw) = -8,(C,0,) 
Ax, (EW) + B5, (C, |) + Те, С) = 5,68) (V-2a) 
DR, Consider, now, that one wants to trace the "constant Xa foci. 


for W =W ith either = É г = TF for a contin- 
n ni i nr δι Un к С) Со) с a 

uous or quantized solution (this makes no difference). That is, the 

following set of equations must be solved, "n" times over the "С" range 


desired: 


for, ow) + вв бо) & Te C90 - -5 €] "7 Š, 


(V-2b) 


каза тала кала тан e 


бо опе must solve "n" times for the intersections of two planes in the 
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ABU space. In other words the constant ШЕ d loci is a set of "п" 


straight lines, with n going toO if Ü is varied continuously in the range 
ο τὸ. 
З. But by letting B = O, for any arbitrary С = fa from the above 


range, in both of equations (V-2b) and solving for: 


η N: αυ к aD 


N 


1) 


дог (С ага) и 7” m 


ENSNoSts a point P, in the OA! plane which belongs also to the desired 


1 


solutions of V-2b for f = É 


ee (see figure 5-4a). 


Analogously for A = O, one сап get a point Р,, on the OPI 
plane, which point will also belong to the desired straight line solution 
of (V-2b) which after all is completely defined. 

4. So this solution can now be projected in the descriptive geo- 
metry reference plane (Figure 5-4b), and mapped by its two straight line 
images. These two images define any point that happens to have a pair 
Ln ríormance characteristics D = a and = ЛЕ ГАК: ' 
for example a point Px in figure 5-4b with A - 1, B = 1, T = 4 would meet 
those requirements. Some work on this has been presented by Thaler and 
Cadena. 

ae The fortunate fact to note though, is that one can keep track 


and interpret the meaning of the two points P P. as part of a simple 


24-22 
patina te solution for шел Ua and ( = oa by simply locating the points 
р ‚ апа P," on the reference plane, since Ро" апа n can be 


located by drawing the normals to ОГ for each of them. Опе should of 


course indicate which pair of ря and. G np was used for this set of 


| 
points. 


6. But now looking back again in the whole process one recognizes 
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that what is necessary to solve this problem is to solve simultaneously 


for each set of values of w =W , and É = Ú a two dimensional, 


arbitrary 
linear parameter plane problem in A,l and B,l and to locate also simul- 
taneously the two corresponding points on the same planar plot. This 
process will continue regularly in order to complete two individual planar 
grids of constant M and constant € curves as shown in Figure 5-5. 

Ue It is important to plot each of the two grids with different 
ls indicators, eg for the OBI and for the OA projections. 

8. If, after the two parameter plane grids are plotted, one wants 
information about the Qa c Uu, É = С), А,В, space locus, the procedure 
15: 

а. First have in mind that the plot obtained lies in the 
reference plane. 
D: Second, locate the two points for each grid, corresponding 
to the E С + Er loci. 
em Draw the two normals on Of to locate the associated 
second projection for each image. 
d. Draw the lines £' and £" which are the images of the de- 
sired locus. 
e. In order to read one of the infinite sets of A,B,l values 
that give шол Б and € = С: 
(1) First, select a Г value, say I - Og 
(2) To get the other two parameter's values, sketch a 
line normal to gO, at the point g. Then: A = ga, B = gb, the intercepts 
of the line bga on the projections £" and £', respectively. 


Dis An example of an actual simple linear-characteristic equation 


system follows. 
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D. EXAMPLE PROBLEM 


8. The closed loop transfer function for the system in figure 5-6 
is: 
i Г! (5-1) » 
во 2 1 1 1 i 
(= + 35 + 2) + Г' (5-1) (А'ѕ + B') 
A СИЕЗ) A  ___... 
СОСО авта) 
where 


n = Sandip = D se l. 


, 


Z: So the characteristic equation for this system being of the form 


previously discussed is: 


SS - s) + B(s - 1) + MS + 3s + 2) + (s - 1) = 0 (у-3) 
Du Practical Considerations on Applying Parameter Plane Methods 

At Chis stage it is ot importance to notice that a digital 
computer program for solving and drawing the constant Да and constant ζ 
loci is absolutely necessary. One could use either a program based on 
the equations of paragraph E of this section (V), or could use twice a 
program designed to solve and draw for the case of two variable para- 
meters, and then superimpose the two drawings (one for an AOl -plane cut 
and the other for a BO -plane cut). But even if one uses the two vari- 
able parameter plane program, difficulties may arise in selecting the 
proper range of A,B,! (that is the scale of the graph) basically for the 
Mésired range of "w" and "f". However, even if one or two curves are 
drawn, the task then is even simpler. 


*2 
To get a first rough estimate of the appropriate range (under 


at. 


4” 


Refer to the conclusions (section V-G) for comments regarding 
root loci methods versus parameter space approach, and also note that 
the solutions by the root loci techniques in subparagraphs "a" and "b" 
to follow, are just particular arbitrary solutions and are not likely to 
be the desired solution. 
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the assumption that some test runs with the computer program of Appendix 
D will follow), one can proceed in many ways. Three of them are as fol- 
lows: 


a. The Case of Forward Loop or Feed-Back Compensation with the 
Forward Gain as One of the Variable Parameters 


(1) Draw the open loop root locus with no compensation 
unit included 

(2) Acknowledging the fact that except for the forward 
gain, the rest of the variable parameters, once fixed, define possible 
compensation poles and zeros, try to find a combination of these para- 
meters that would change the open loop root locus to be stable for the 
desired working frequency with a Ç approximately equal to 0.707. 

(3) As a very quick check on the approximate forward 
gain range apply the Routh-Hurwitz stability criterion for the closed 
loop characteristic equation. Sometimes this is not enough and one should 
draw the root locus. 

(4) The rough scale to be used for the first time should 
be as a rule of ος ten to a hundred times the values of the com- 


pensation variables selected and the maximum stable forward gain. 


Da Compensation (Forward Loop or Feed-Back) with a Specified 
Forward Gain 


(1) Insert a variable gain parameter in place of the 
specified forward gain value. 
(2) Repeat the procedure for the previous case, but ac- 


cepting the resulting scale if the specified gain is as a rule of thumb 


d 

53 Tne rule of thumb proposed, might of course not. work 2. 
very large or very small, and therefore some change should then be 
necessary depending of the particular system. 
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twenty percent less than the maximum stable gain obtained. Otherwise, 
try a new conbination of the compensation variables. The reason of this 
procedure is to assure an A,B pair which would give under the operating 
Mat least one complex pole pair in the region of 2 = (.3 + .8), Ъе- 
cause the possibility of not getting a constant € curve is reduced when 
MS palr is forced in the corresponding s-domain area, for "Àj" around 
707. 

oF The Closed Loop Characteristic Equation Technique 

4 When only the closed loop characteristic equation is 
given, and one has no idea about the actual system, the following pro- 
cedure could be followed: 

(1) Obtain three Root-Locus Plots, each for one of the 
three parameters, as the only variable and the remaining two set equal 
Bora convenient constant (say unity). 

C(om Tror a certain range of "C" and р desired, define 
the corresponding range of each of the variables in each of the root- 
Pocus plats. This will give an indication of the scaling of the problem 
to start with. 

One should also have in mind, once engaged with a time 
delay problem, to set the time decay constant equal to zero, to apply the 
method discussed above. 

4. Practical Considerations on Applying Parameter Methods of 
Paragraph C 
Considering equation (V-3), and having in mind the discussion of 
paragraph C, one comes to the conclusion that if a BOI cut of the para- 
meter space is tried for, the particular system of figure 5-6 (i.e. A=0) 
ms physically restricted from having any С, since the roots of the 


GH function will all lie on the O axis in an "s-plane" representation. 
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So there is a physical limitation in applying the technique, 
which one can overcome by a change of variable А" - A 4 l, thus main- 
taining the tachometer feedback of the system for the case of A = 0. In 


this case, equation (V-3) becomes: 
2 2 2 
Als - s) +B(s -1) +I(s + 35 + 2) + (5 -1) = 0 (У-4) 


9 Having in mind all the practical considerations above, one Can 
proceed to specify an approximate set of A, B, l range for a first try 
of the parameter space solution. 

a. From figure 5-7 (following the procedure section V-D3), 


опе can see a good value for the compensation zero to be at O = 2,(or 


Q 
l 


ου ο). This can result for A! = 2, B' = -4 and equivalently for 


> 
i 
E 
ο. 
II 
i 
л 





b. Since the Routh Hurwitz criterion is going to give an 
expression ОЕ Г! "τν, τυ; rPESRoLbows "that: 
шах gain gain 
l; — ) 5. 
тах Γι gain 
ma x 


mer stability, and this is not enough. (Root-locus is needed). 
Cy For A =.1, B = -5, the second order system of equation 


(V-4), becomes: 


22 + ος 25 + 2 


Г+2 
EL follows then that for al = l//2 and v = Оа а Г = 10 is needed. 
d. So as a first try the following ranges should be tried 
out: 
(1) |А = 10 
(2) |в| < 50 


(3) |T| 5 20 
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6. In figure 5-8, the results of the parameter space analysis are 
shown. In order to show that the results are valid, select the parameter 
space loci for a = 1.0, © = 0.7), represented by its two images (the 
dash-dotted lines). An arbitrary set of values for this loci is, 

(2-059075 952 32 s 12.0); 
which gives a second order polynomial, (equation V-4): 


πι. σος 
ανα ο μμ 5 - 





5 


with: Ше 2955 
' n 
É = του 
For a second demonstration consider the loci for (w = 1.6, É = .6). An 


EuPraryv set of values for A, B, T on this loci is, 
Cee ee ='-).5 T = 3.0). 
giving a second order polynomial,(equation V-4): 


2 o. 215 8.500 





on 
with: See 2 
n 
£ = .603 


These results turn out to be fairly accurate and the theory presented 


OrKs in practice. 


E. SOLUTION TO THE GENERAL PROBLEM 


E; Consider again equation (V-2). The discussion in section 
(V-B) indicated the necessity с an infinite number of cuts to be done 
for the complete parameter plane solution in three dimensions. However, 
depending on the desired degree of accuracy a finite number of cuts will 
Sufsice for practical purposes. 

7 The solution for a case with more than two cuts will be demon- 


rated in section 5. 
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2m The equations for the general cut through the T-axis will be 
derived by letting: 
А = АВ (V-5) 
as the plane of the cut implies. 
m By substituting equation (V-5) into equation (V-1), the fol- 


lowing set of two equations results: 


2 2 S 
BAY, + δι) + Гє, zu; Ad, m Bl (^n, + 91) + В ГАн, 5 (V-5a) 
2 2 re. 
Βίλα, Б 5.) + Ге, + В Ad, + Bl (m, +.) + B Au, Zu (V- 5b) 
Therefore: 
2 
ES "αι + 9,)B AG, B 
= πι. (У-ба) 
ει +В (Ату + 94) +В Ma y 
and: 
4. 2 За | 
BAA, +B SA (A. + AÁ +À( A, +A | + 
LD L As no? ( u Ô E 
в? | ZA +À(A +A, +A А) +A Ja 
тр pa т єў ш δ 
[ | ін С 
В pu + Ang? m (А, ΕΤ: Δε 0 (V-6b) 
where in general: 
х Di 
c 1 1 = r -- 
um r Riu 5199 
D 
5. So in order to define a certain desired ше с point in the 


А,В,Г space, one should follow steps as follows: 

а: Define the coefficients 94, Bjo sees 51 and oe 855. «55, 
which become constants for W = Wis б = С: It is assumed of course that 
the time delay in the system is fixed. 


b: Solve for B from the equation (V-6b). A digital computer 


pug 
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routine will be necessary for this purpose. However, the conventional 
solution of a quartic equation is listed as Appendix E, for easier ref- 
erence. Care should be taken to exclude any complex roots. 

m Then T will be given by using (V-6a), for each admissible 
E root., 

di. ne Ck Ol the cut=plane is kept the other projection is 
not necessary, since А = АВ (V-5), and therefore for each point in the 
ВГ grid that will be generated by the procedure of subparagraphs 58 
5c, there will correspond another one as shown in figure 5-9. 

Im figure 5-10a the grids of constant É and 4, curves 
obtained for three different cuts are shown. Notice that one projection 
K васп спі is sufficient. 

To construct an approximate constant шие С. согхез ове 


Mist select the points РІ, Р РА on the corresponding grids ( one 


2» 
from each available grid), and by joining them by straight line segments 
obtains one of the two required approximate projections. 

In figure 10-b this projection is shown, with the con- 
Struction grid eliminated. ‘The second projection of the constant Ю 


5 curve can be approximately defined, by locating the projections 


А À ete. 


Ж ҮЙРЕ Pl from equation (V-5), with A = A ›› ^з› 


27 3 
Е. EXAMPLE ON A SIMPLE DESIGN PROBLEM 
Consider the hypothetical system of figure 5-11, and suppose that 
the operational requirements of the system are met if the operating point 
is W = .4 rad/sec and ( = .5. 
α Q 
Consider also that the cost for building this system was found to 
ET 
be C - M: where M is a big number and T is so small that e 5 = 1 


(in order to simplify the algebra). There exists no information about 


DEE 
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Ehe Operating cost of the system as a function of A, B, Г, зо this aspect 
will not be considered, but of course the stability is presupposed. 

ШЕ 15 desired to select А, B, T that will simultaneously give s 
E. as specified and minimum possible cost. 

It is noted that A physically represents the positive gain of 'an 
amplifier, therefore there exists the restriction A > O, 

SOLUTION 
E. The closed loop gain of the system turns out to be: 
Ts 


E LT (s+3) (As + Ше. 
Р (54% (2.828 NS (442.8288) 5^4 ABs | 1e 75 Ar (s^4 3s) + Г (5+3) | 


СУ 


So the characteristic equation for the general cut-plane A = AB, and 


ШЕЛ will become: 


νη] 
Seis) = D + 2.8288” + 45 | ο. + 2.8285” + s | 11] 539] + 


L 


2 i 
1-55] - ο. 
L. - 
2. At this point one must make а decision on selecting the number 


 СПЕв”Сһасгс will suffice the specific application. For this case con- 


Sider eight cuts, corresponding to: 


ЕО Е τι uo Neto 
In table 5-1 the resulting CE(s) are tabulated. 
E By applying the method of section VE, the curve for constant 
(ο, = .4 and € = .5), was obtained and displayed in figures 5-12a, 5-12b. 
These figures have been plotted according to the data of table 5-2. 
In figure 5-12a the interesting situation for the stable com- 
pensation of the system under design 15 displayed. 
4. According to the available information: 
a. Gain Upper Limit: 1=0,217 units (stable operation on 


specified Ws C). 
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b. Gain Borer Шш pp ος. Nevertue less rl cas possible 
to. get a minimum gain of .184 units for the lowest B parameter employed. 
ей Compensation Pole: 
0.00284 < В< 0.47 in the range Í = (0.184 OZ 
die Tachometer Feed-Back: 
090050755 A ~ 2.84 in the same | range. So it follows 
τ ἔπε system is more sensitive in "B" than іп "A" variations (for this 
specific I range). 

5% So in order to satisfy the requirement of as small gain as 
possible, imposed by the specifications and having the information of 
figure 5-12a, one should choose of course the set: 

(A, B, Г), = (2.84, .00284, .184). 
But if it was possible to pay 9 percent oped the solution 
(АУБ; 1), жЕ о 220) 
is much more advantageous and takes care of possible component drift. 
6. n ubscstutipnge the set (A). B, |) € (6.115. 25. Ὁ in CEs), 


(for T-0), one gets: 


я 2 ` Г 2 С 
Cis) = B5 + ο... + 4,707 s + 1.005 | + 1.155 + 4.445843 | = 0 
the root-locus plot of which gives GM CU еее 56, ΠΠ του 
Í = ,2, which is very close to the desired specifications (see figure 
5-12с). 
l: It is of interest to note what the constant G = O and € = 0 


surfaces look like. 
а. From CE(O) for s = (σε 0) the O = O surface is the 


Ш= © plane. 


3 А Г -Г 
ES πο е не cost function is C = M it follows that =. 5_—1..087 


Ы. 
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b. From CE(jw) for the s = jw surface one gets the following 


set of equations. 


ЕГЕ ρω - AFu2 + 3P = 0 


- (2.828 + В)? + ABO + 3AF0 = 0 
ο By simultaneously solving the G = O and € = 0 equations 


Bezequation for the cut of those two surfaces is: 


p^ + 2.828B + & = 0 


which gives: 
Ве - /2 4 j/2 


Bherefore there exists no intersection of the Oo = 0 , Ë = 0 surfaces. 


em CONCLUSIONS 

From the theory presented in this section theoretically and shown 
to actually work for the simple system of figure 5-11, the following con- 
clusions can be obtained: 

15 Ше ив possible to constmuet the constant "W-C" curves for a 
three variable parameter system. Care must be taken in order not. to 
exceed the physical limitations of the system, e.g. taking cuts that 
would perhaps destroy the feed-back om forward gain of the system. 

2% In the a the variable parameters show in a nonlinear 
fashion in the system, the constant 'w-£'" curves have no straight line 
Шиодесетопс. 

з, Care should be taken to differentiate between different seg- 
ments of the solution. It is of importance to define which of them 
provide stable and which unstable operation. Until a well grounded 
stability theorem is developed in the three dimensional space, this task 


which was easy for the system of section (V-FYmight prove impossible 
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for a complicated situation. In this case the method of this section 
will be handicapped. 

4. Another thing to notice, is that for the nonlinear са$е, (А,В,Г 
Products), a is possible to come up with a curve which is more or less 
normal to the l-axis. If this happens it is not easy to read the appro- 
priate A,B,l sets in one hand, but sensitivity on the lI parameter is 
proved. So if one wants to operate under the conditions imposed by 
the curve, one must set I = Ци and then possibly one can reduce the prob- 
lem to a two dimensional one. 

В, Sensitivity analysis of a system is possible for a multivariable 
control--system, by applying methods described in this section. This is 
in contrast to the approach of using a family of root loci where the 
method of sensitivity analysis is not so direct. 

6. Another great advantage of the parameter space approach over the 
root loci method is that system optimization through adjustment of para- 
meters is more direct, and illustrative. 

7. Of course in order to arrive at the desired result a vast 
amount of analysis was necessary, and it is a disadvantage which of 
course is removed when a computer program similar to the one in Appendix 


Dis available. So indeed, the method is "digital-computer" oriented. 
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VI. AN ACTUAL ENGINEERING APPLICATION 


А. GENERAL 


E In this section the roll angle control system for a vertical 
take-off aircraft is studied. 

2. in 2 post-design analysis of this system, (which was initially 
designed by the root-locus methods and thereafter operated successfully), 
attempts to locate the actual operating point inside the two variable 
parameter stability region were originally unsuccessful. This implied 
that the stability criterion of chapter IV was invalid, and since the 
derivations of the stability criterion was rigorous, further study of 
the particular case is needed. 

9, In the subsections to follow, the system is first considered 
from an elementary aspect and since it needs compensation, its possible 
modes are defined. 

Then a parameter plane design follows and having in mind the 
pesults of the stability analysis, an appropriate operating point is se- 
lected. Time domain simulations follow at last to insure the correctness 


of the parameter plane results. 


В. THE SYSTEM 


1. Brief Physical Description 


а. The system is presented in block diagram form in figure 


In the component analysis of this figure one can observe 
seven variable parameters. But since the servo constants have been pre- 
set to their values of table 6-1, the only variable parameters remaining 


are K , K , and T(which is bounded within limits). 
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TABLE 6-1 


ΤΙ OF CONSTANTS USED IN THE SYSTEM OF FIGURE 6-1 


CONSTANT VALUE 


T ‚15 + „30 

K 1/8.1 = .1235 
D EO IS 

E EUM OPE 

F 2.329 
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p Since the permissible range of variation on T is small, it 
ENrcasonable to start with the minimum expected value of T, (0.15), since 
the results of section IV show that this should give the maximum stability 


region in the A,B plane. 


2 Ob jectives 


a. To actually define appropriate z K, values for stable 


ф 


operation. In doing so one should have in mind that zetas around 0.7 


and frequencies from 1 to 3 rad/sec are desirable from a "pilot-perform- 


*5 
ance" point of view. 


7 


Bis To check and find if the operating values of Е K, lie 


Ф 


in the parameter plane stability region in order to check on the dis- 
crepancy reported (section VI A). These улпс калге: 


) = (95, 190) 


(К Ka 


з, Sestems Transfer Characteristic 


Ты. 72 


к. | 7 
/ 
$6, = (.305) De “is + .444s + 488) ED 


F(s) + с779 дер, (е) 4 вер, (6) | 


бы 
| 
E 


where: 
722 2 
А“Р, (8) = A. s + .666s + .4865 + 0.0] 
НЕ - 
us 3 2 1 
В.Рр (5) = В. 05 +5 + .4445 + .488 | 
and: 
T (2305)K K 


=a K > Tachometer 
D р р 


В = EUR K. = 77.2 = = Proportional Control 


D ф 
6 5 & 3 
F(s) = s + 98.46968 + 4941.16808 + 8978.0825 + 


5086.9961s* + 3162.60115 - 147.59 


m 


^ 


These were established by the original system designer. 
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* 


a. The Systems Characteristic Equation is: 
6 5 4 3 
CE(s) =|s + 98.4696s” + 4941.1680s + 8978.0828 + 


5086.99615 + 3162.60115 - 147.59 | T 


ΛΑ ES B Es + .488s + 0.0| + | 
a B | 0-65 -- я. + .4445 + 488 | 


= 0 (VI-3) 


с. POSSIBLE COMPENSATION MODES 
Ша By looking closer to the feed-back component of the figure 6-1 


System, we can see in general that its effect on the open loop transfer 


` 


function is the insertion of a compensation zero "2" 


паша by controlling 
the values of the variable parameters iS and Ss one can cause various 
compensation modes. 

2: Noting also the results of section II, where it was observed 
that though the shape of a root-locus was symmetrically deformed towards 
ΙΙ», its basic form remained unchanged with varying delay, one 
can search for the modes mentioned above using the delay-free version of 
figure 6-1 without loss of generality. 

Sr In figure 6-3 the open-loop root locations of the system are 
shown, with compensation 22 not included. 


So with no compensation ος ρουμ ο produce instability in 


1 


the desired operating range. 


Море абсо? first, that poles Ες апа Ες located at a natural 


frequency contour of 69 rads/sec with a = .7, are too far away to 
be considered at the desired frequency range and second, that zeroes 


2; 2, restrict root-travel branches leaving poles P to remain 


~ 


in a very limited range in the s-plane, so that they force the continuous 


er 


presence ofa complex root pair with a set of (£, X С 7 окпе 


average. Theoreiical cancellation would of course happen for infinite 
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gain, but since these four roots are very close together in pole-zero 
pairs, the influence of (P4, P,) in the systems transient characteristic 
will be much reduced. 

4. | Now the real zero p) introduced by the feed-back, can general- 


ly take three possible relative locations with respect to the P ΕΡΕ; 


Шон? 
as follows: 
а. zZ ecorthe Lett ’oi P_. 
С 2 
ШО сыг eause a Closed form of the root locus and for a 
MProper gain will provide a pair of complex roots in the L.H.P. solely. 


Ф 
b. 2 to the Right of P 


In this case, both is and K, must be positive. 
Е 
In this case one of the two feed-back gains must be nega- 
mye. oO if in the feed-back expression 
Bor + Kg? K, is factored out of the parentheses 
1 
опе will have: 
K 
K (s E 
p К 
Р 
which implies a negative gain locus, for га < rO; 
C Z. Between Р, апа Р,. 
ШЕ ο οἱ теси ро a гез! root pair, one of which will 
be in the right half plane, for "low-enough" gain. 

Die Mer this discussion,it is) obvious that two are the possible 
compensation modes, one with the ne Iingtehes Her. and one with 422% іп 
the R.H.P. Note that the operating point of the system was regulated for 
a Left Half Plane compensation with: 

(K , K) = (95, 190) 
p 9 


 РІЛІпе а zero at е -2. 
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6. In figure 6-4a, the root locus of the system is plotted for 
this selection of SE DE This was done as usual by equating the open- 


loop transfer function to minus one, and thus getting: 


6 
Е + r. + 5056.38" + 8966.6053 + j + 


' E Е 
4905.43s - 225.3 | + ο] ο E >. + 1.371s + :973 | = 0 
where G was the variable gain parameter: 


і 


- | = 59400 К, 


р 


ο. К x K = 
Eel NES Ыра: 


Eu 


and since the working K of the system was selected to be .1235, one gets 


the operating point to be as follows: 


G K K K А P С UN 

P 2 rad/sec 

Е 320 61235 08 190 7340 14680 „812 Ше бо 
T: To illustrate by the root-locus technique that compensation is 


possible for 2 in the Right-Half Plane, one can select Sr Dann 


and get from the open loop transfer function: 
Б > ОО ΠΠ υπ. 


4905.43s - 225.3 | ΠΗ - Ee + .O44s = 497] = 0 


where: 
πι. X K 
Els Р | = -59400 К 
L о (K, = -95) 


So for positive K, as in this cass a negative gain locus is resulting 
and is displayed in figure 6-4b. In this case, ee for stability, is 
equal to .0965 (equivalent for G' = -5722), which is a restriction for 


the designer. Notice also that for "{" in the area of .707 or greater, 
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the natural frequency will be less than .55rad/sec. 
Eo. With this information about the system, its preliminary investi- 


gation is concluded and the parameter plane design follows. 


D. DESIGN IN THE PARAMETER PLANE 

1. Ву using the computer program of Appendix D, and the selected 
sets of constant "E" and constant "ш" curves shown in table 6-2, three 
sets of parameter plane grids were plotted for values of the time delay: 


Оа ОБИ э 


I 7027 
These results have been plotted in figures 6-5a, 6-5b, and 6-5c. Quasi- 
linear versions of the maximum root regions found T chapter IV have been 
superimposed for easier E EAN 

2. Having in mind the general range of (С, ) desired, (section 
ШЕ) апа the particular (C8 0) = (.812, 1.65) that resulted to be the 
ones preferably selected by the designers of the system, one can plot 
the variation of the (ο, = i O and 7.8) Constant (w »Ὁ) curve for 
varying time delay. This plot and its expanded scale version both re- 
sulting from data of table 6-3 are shown in figures 6-6 and 6-6a. Notice 
Bie approximately linear relation in these plots, especially within the 
operating time-delay constants T E [.15, .3] where AT and displacement 


in the A,B plane remain in a constant ratio. 
TABLE 6-2 


Constant "Ё" 
Curve Values 


GE ο ο ο ai C= в С. - 85 


Constant "w" ¿E TA _ 
un c ues E аа w, = .688 Ww, =1.0 ш, =1.3 ως Ξ] .6 


This fact can help in predicting appropriate A,B values for opera- 


tion with a "time-delay constant" different than (T, T5, T = 2 gm 


3) 
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However since there exists no limitation on the lower value 
of the time-delay to be used, зоне other words since ө = .15 is an 
admissible delay constant, the best choice of A,B values, in the absence 
of any other criterion to be met, is the set for T = ο 

ПАУЫЛ - (25742 9453) 
B. After having the first complex pole pair fixed, one should 

think about the other two pairs. 

а. Іп the case the system is operating within the stable 
region, one should take into account that the complex pair at esie 69 
rad/sec contributes a very small percentage to the total system transient 
response and therefore can be neglected. 

Ь. 1f the second pair remains near the (0, у= (у 55) 
zero pair, in the case now that finite delay is present, there will again 
occur small interference on the transient response from P,, Pye But in 


order for this to happen one should be able to observe constant (x) 


curves of about (.7) all over the stable region for time-delay constant 


nenes range from (.15 = .3). To show that this is true the parameter 
plane curves were solved for an average time delay value, (T = .225), and 


plotted in figure 6-7. So by observing this figure one can see the ex- 
istence of the D 0.7 curve spreading out across the A,B plane and 
therefore present for every (A,B) selection. 
4, As far as the solution which was obtained by conventional 

means (root-locus, etc.) one can observe in figure 6-6a that, for T = 0.0, 
the set 

(A,B) = (7340, 14680) 
corresponding to 

e К.) = (95, 190), 


Ф 


lies in the statle region. 
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VARIATION OF "A' AND "В" VERSUS "T" FOR CONSTANT OPERATION 


@ 15 


25 


units 


7340 
2098 
5722 
5574 
5280 
4893 
4463.6 


4026.4 


TABLE 6-3 


ON (0 6) ООВ) 


units 


14680 
2062 

10568 
9453 
8279.4 
730 
6064.1 


511023 
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E. TIME DOMAIN SIMULATION 

T. It would be convenient to present at this stage an appropriate 
transformation of the system from its physical block-diagram to another 
mathematically equivalent diagram which is appropriate for simulation. 
Two such mathematical block diagrams are shown in figures 6-2a and 6-2b. 
The conventions used in figure 6-2a are shown in table 6-4. 

2, From this time-domain simulation the time responses of figures 
esa, to 6-3f were obtained, with the results as shown tabulated in table 
6-5. 

3. The steady state output one gets from equation (VI-1) is given 


by the formula 


ñ 36.7 
ЖЕ EN ορ шу 


Б, CONCLUSIONS 

L: NUGESvSten proved to be stable for T 9 .l, .15, .2, .25, .35. 

D. The results of section IV and the time simulation coincide as 
ШІГ ас (іс Т - .1, T = .2 cases are concerned. 

БЕ Kor che DT = .3 stability analysis, it seems that for both the 
section's IV method and the time domain simulation, one would expect to 


get stability for T = .3 by trying bigger frequency and final time ranges 


respectively. But since it was decided to use T = .15, this was not 
necessary. 
4, The system is safe stability-wise to use for the selected 


parameters of: 


А = 5574 = = 72.2 
B = 9453 E - 122.5 
with a time-delay constant of T = .15. 
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The time- domain characteristics for a step response are for 


пиле case: 


" Eu OLE DUE - 19 
MeS O y State "Output" i 


Ew πα for "Max Output" = 9,19 sec 


Y. Rise Time (107-907) 4.86 sec 
0. Time to reach (1%) error= 10.4 sec 


eee Error at "t=20 sec" = 2.8 ppt 
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TABLE 6-4 


le. Delay 


- Ts 
e = T (s) 


D. Numerator 


Gall I 


100 





CN п 


17 (.488) = L, 


=, 
ll 


Lx T (s) 


z 
lI 


1 Ly кх T (s) 


N. = L, ок T (s) 


B. Denominator 


1.0 


а 
Il 


D. = 98.469 

D, = 4941.1680 

De = 8978.082 + A x T (s) 

D, = 5086.9961 + L(.444) x À + В) хт (в) 

D, = 3162.6011 + Г(.488) ea) x Bl x T (s) 
D + -147.59 + | (.488) хв| x T (s) 


0978082 


O 
O 
Lo 
II 


DD2 5086.9961 


DD1 3162.6011 


1] 


DDO = -147.59 
4. With the conventions of table 6-2 at hand, the system's transfer 


function becomes: 


E _ N. (s)s + N, (s) 5! + N (s) | 
re REC E 5 4 7774 2 | 
6 s)s + „(s)s + D, (s)s T D4(s)s + D, (s)s + D, (s)s + D (е) 
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or 


1 1 1 
> ms S EG i No 6 
EN s S S 
R 1 І 1 1 1 l 
l + D. ЖЫ: D, 5 + b, ER р, τα D; gt D. EG 
S S S S S 


From this form the signal-flow diagram of figure 6-2a, follows easily. 
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APPENDIX А 


ADDITIONAL INFORMATION CONCERNING THE ROOT-LOCUS ANALYSIS 


OESCACLOTSOTBHOWITHOCELHE-DEBAY 


This appendix contains more specific information concerning the 
7 . . . . : 
Pade approximations as well as data for the characteristic equation of 


the system: 


t 


t 
Т" B (see Table II) 


resulting from the block diagram of figure 2-2, but with the above 
mentioned approximation inserted. 

The Padé Approximation 

If one represents by DE(s), the transfer function of a time-delay 


in the "s-domain'" then for T being the time-delay constant: 


2 3 
DE(s) = e τοσα τη Q2. -AR πι 


As already mentioned in section II, the objective of the Pade ap- 
proximations is to define two polynomials N, „(Ts) and D. „(Ts) of order 
2 2 


να and "b", respectively, which have a quotient 


N (Ts) 
F p (Ts) = а 


ум 


such that a certain number of terms of both the related series expansion 


Of DE(s) and this quotient are identical. 





Letting: 
N_ (Ts) A, +A, (Ts) +) Е πη 
- š О 1 2 а 
ЖК е 2 b 
4 a,b Во + u, (Is) +4,(Ts) EH y (Ts) 
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one observes that ВЕ y (Ts), icta function ol rhezordorpa b of No y CIS) 
ә 3 
and D. p (Ts). So after selecting a,b, one must choose the appropriate 
2 


values for the coefficients ic = ^ and ТИ = μι, in order to define the 


£ - : А 
Pade approximation. 


The polynomials N. p (Ts) and D, y CI) are specified by: 
` 5 


3 


[=a m m=n- 1 
x (Ts) a-m 
= iS A AS 
5) ; у 1) ща | a+b-m 
nel m=0 
п=Ь n ΠΞ -l 
E alles) П b-m 
D ^ ТД г! а + Ὁ - πι 
nsl πο 


In Table Al, a list of the first six Padé approximations is tabulated. 


Ше 





LIST OF 


ШАП! Ts 


TABLE Al 


THE FIRST SIX NORMAL (a=b) AND REDUCED (a=b-1) APPROXIMATIONS 


Р, 109 


2 


Е, 20) 


2 


Е 229 - 


3 


K ко) 


Note: 


Б 
ee x 1 

- + x ο л 
12 - 6X 4 X* 6 - 2x 

F 2 F, 500 = 2 
Dux : ООУ 

| 04 2 

ОООО „12у l X O 

E 2.3 Fo 36%) = 2 3 
120 F OOX + 12X + X 7 60 + 36X + 9X + X 
1680 - 840X + 180x* - 20х? + ха 


1680 + 840Х + 180x* + 20 + x" 


840 - 360K + ох“ Е AX. 
840 + 480X + 120x* + 16x” + x" 


30240 - 15120X + 3360х2 - 4207 + 30x" - 
7 
ο ο 51T F 3360X 420% + 30x* + x 


Шота 5720 4 1260Χ΄ - 1200 a 5x! 


15120 + 8400X + 2100Χ΄ a 300x> κ 25x" + x? 


665280 332640Х + 75600x° Е 10080x” T 840x" - 42x? + x° 


665280 + 332640X + 75600Χ΄ ЕН 1008ох” + 840x" + 42Х? zs x? 


332640 151200Х + 30240x" = 3360X> Sr 210x" - өх? 


332640 + 181440X + 45360X^ + 6720Х? + 630x" zs one x 


6 


For more detail on the Padé approximation, see reference [3]. 
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TABLE А2 


DATA RELEVANT TO THE CHARACTERISTIC EQUATION OF THE SYSTEM SHOWN IN FIGURE 


ο. SECTION ΤΙ 


Bharacteristic Equation: 


r / 
(0)89 + (1680)s- + 1840Т 2 + 1801 А + 2017 EE T > + 6 + 
i Е И um n 4 
1680 840T 180T 201 Т 


l 
© 


K | (1680) 89 - (840D)s X (18072)s2 - (20T)2)s 4 (ye | 


Galeulations For Repetitive Terms: 
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APPENDIX Б 


РАКАМЕТЕК PLANE (CONSTANT RM and ''C'') EQUATIONS FOR A SYSTEM WITH A 


TIME DELAY 


чм 


The characteristic equation of a system is given to be: 


n X 
ИС з κ. _ 
pu S + а с 5 Ξ 0 (В 1) 
k=0 k=0 
ШЕ: E 
where ар АБ + Bc. (B-2) 


These are the equations (III-2) and (III-3) which characterize the 
system of figure 3-1 encountered in chapter 3. 


The objective is to derive a set of two cquations of the form: 


в] 0) + BC (20) δις) ο 


A*B,(9 ,C) * B*C, (o ,C) + РФ) + 0 (B-3) 
r = 
where: ШЕ. с; р, | БОК аср сг ГС Оп ot Y... and C 
ша 1 ij. n 
i=1,2 
and A,B the unknowns. 
Using the fact that 
= (cos Ө > Sin 9) =W К. 
EAG j zZ 
it follows that: 
4 | 4 x 
IK As ein) 39. - 
ае $ = /ау е п we ) = 
k=0 k=0 
А. 
Уан е7 ТХ Cos 8 - jT sin Ө 5 K „JK 2 
К п 
k=0 
£ 
B K -Tuncos 8 j[K9 - Tu, sin 0] 
y at n n 
арш е е 
к=0 
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jl ko - Tu sind] “τη σος ο 


Now expanding e 


and getting e out of the 


summation Operator, one can get: 
£ | 
m -T 7 
Зан е ES (B-4) 
k=0 


w | cos (Ke AS A SA Ty sin 33 


Emularly then: 


n n 
NU INN к 

е ^" : STR _ 
]94 5 24,9. (совКӨ -- jsin K9) (B=5) 
k=0 k=0 


By substituting equations (B-4), (B-5) into equation (B-1) and by separa- 


ting the equation into two parts, one for reals and one for imaginaries, 


` 


ШЕ то Томс: 


п L 

N K -Tu_cos 8 5” K 

y е > n у " 9 c n ° = 
ie, cos КО + е DEM cos (KO Tt sin Ө) 0 
κο. 0 

п $ 

eei y -TE ха = 

/4 МЕ oda пещи ан МЕ (е en Bye 

Коп n n 

k=0 0 


and by substituting 


ар = AD, + Ρε.» (6-2), one arrives to the following two equations: 


£ 
-TUncos Ө T Т K = Ч > | 
ale Rn соз (КВ Tu sin 6) + 
к= 
$ 
-Tu cos ð` K | | 
Y = 
B [ο n erts cos (КӘ Tu sin Ө) E 
К=0 
п 
| /4 Ὁ К COS кө | η (В-ба) 
K п 
k=0 
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and: 


ү 
ШОЛ сов б s K .. - | 
А | е pe, sin(K9 - Ty sin 8)| + 
k=0 
4 
EE ο cos ON К 1 
( i - Ty = 
B E Se, ѕіп (КӨ То sin Ө) | 
k=0 
n 
í Md M ] 26 
| Де sin КӘ | - 0 (B-6b) 
k=0 


One should notice that equations (B-6) are more or less of the 
form of equations (B-3), which are the target, since cosé, sind, апа Ө 
can be expressed as functions of €, so theoretically, the problem is 
solved. But since in general a digital computer solution is desirable, 
it is noticed that a summation in K containing 

sin(K9 - Tx sin 8) or cos(K9 - Tr sin 8) 
will not be too easy to evaluate. 

So this is the point that one must refer to reference [4] ana 
appreciate the fact that he found some existing recursive relations in 
his parameter plane analysis without a time-delay. From this fact, ES 
of value to think that perhaps the same thing applies for the case under 
consideration 

For this reason one can try: 


le Firs[ Set of Substitutions. 


Call: 
Yk SUN = Tx sin 8) (B-7a) 
X. = cos(K9 - Tu sin Ө) (В-7Ь) 


ме 


From equation (I-7b), by substituting соз Ө + -С, sin O Eze), and 


by appropriate trigonometric expansion it follows that: 
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24 
ll 


ar ары, cm [o 2 


PS 
N 


0 cos Са 


1 πα sin E ο. | - b cos Im Ae | (B-8) 


< 
I 


Y V pu (His iS nort 


Similar relations can be derived for Y i о 1 


necessary. 
ТЕ can be shown by induction applied to equations (B-7), (see Appen- 


Huc). that: 


ic οι (5-28) 
. "Y. i E (p-9b) 
1 » 
Y. - m LEX. 4 en (B-9c) 
s 
2 Second Set of Substitutions. 

Call: 

U. NES KO (В-10а) 

Т = сов КӨ. | (B-10b) 
ШООГО Томс that: 

ΤΙ ὁ 

10 = 1 

τι = -C (B-11) 


and again by induction from equations (I-10) it can be proved (see Appen- 


uc). that : 


Te 7 "231. - T (B-12a) 
L ---. κι (в-12Ь) 


1:57 








ү! 
Now after noting that 
Een cos 9 _ E 


and by using equations (B-7) and (B-10) into equations (B-6), one can 


get equations of the form (B-3) with: 
SÓ 


σι = CQ ы X.) 


D, = D (w> ар T.) 


1 1 
τ" 
С, = Cr, С, Y.) 
η (В-13) 


ПЕСИ па step in this derivation is to use equations (B-9c) and (B-12c) 


which provide Y. and Ur in terms (Xy Х 1) and CE le] 


E nally from (B-6), (B-13), (В-9с), (В-12с) it follows: 


) respectively. 


2 5 
А ет һы Kx qe, uj +в [ета бё n | 
E 
12 де, ΓΠΣ, | = 0 (В-14а) 
4 
А p лыш Қыш) d XC.) | ah 
K=0 


в (е ΠΟΤΟ 
К-0 


| ао (emo t 0)] = 0 (B-14b) 
K=0 
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where: 


ο πο ο (B-12a) 
Mea 7 2% eat (1-93) 


Equation (B-8) and (B-11), giving X x ande 2, T are 


ο i ο 
also necessary if a computer solution of equations (B-14) is implemented. 


The definitions of LB, С.» Dilii 2 follow by comparison of equa- 


tions (B-3) and (B-14): 


3, 
' U „Tu 2 = x | 
В, = JA a W “x, CW) 
K=0 
5 
е y 
E Er ο, RC, 
К-0 
п 
0, = day j| TG) 
ENS 
4 
Ж Шер“. т К 
20 s PO LEE) + Хе р, wed 
К-0 


4, 
Co n enS еа сха) = 5.1» UE 


K-0 
n 
D, - dr ltr C) * T, 0] (B-15) 
K-0 
К 2 З SE 
БО IÍ BC, ВС, 0 one can get from (А-3): 
τυ, “υπ, BC 
в = (B,D, - D,B,)/(B,C, - B,C,) (B-16) 
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APPENDIX C 


PROOF OF THE VALIDITY OF THE RECURSION RELATIONS USED IN APPENDIX A 


г. To prove that the relations (I-9) are true one must first 


prove their validity for K = 1, and then given the truth of these rela- 


tions for K= V one must prove their truth for K = УН. 


p validity o  Relations' (B-9) for K = 1. 


а, In that case it is required to prove that: 
(1) X, = -25X, - Χρ 
(2) Y, = -2бҮ, - Yo (C-2) 
τ. 
ME) no (C-3) 


mayen, Of course that, 


K 


and X, 


Б, By direct substitution of X 


sin(KB - Tx sin 0) 


cos(K8 - Tx sin θ) 


Xo X and appropriate 


ZEN 


trigonometric expansions in equation (C-1) one gets: 
cos 20 cos (Tw sind) + sin 20 sin(Tw sin Ө) = 


(-26cos0 - 1) cos (Tu sin9) - 26 818 sin(Tu sin Ө) 
and taking into account that: 
G = -cos9 


Ms Obvious that: 


Hi 


-2% cos ð -1l = 2 сов т = σσ. 


sin 2 


ІП 


апа -2C sin 9 = 2 sin 9 cos Ө 


Therefore, equation (C-1) is true. 
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cn In a similar way one can get that equation (C-2) is true. 


d; Substituting for Y X. in equation (C-3), one gets the 


ΙΙ 


following: 
sin 8 cos (Tu sin9) = cos 9 511 (То 5119) = 
pe 
(1-С) ‘LEcos8 *cos (Tx sind) + Csind "51 п (То sind )+cos(-Tx sind )] 


This relation is true since € = -cos 0 as mentioned before. So equation 


509) is valid. 


B. Now in the case that K=V, one must prove that: 
ο ο ο ES (6-7) 
(2) Yoyo = EY - Υν τ 
o 
e T d O) [ex + xO] (С-9) 


provided that 


Rum XV ул (0-4) 
τ... en 
, κ ;.- Πο) πο. r И (С-6) 

4, The proof for the three sets of equations in pairs of two 


follows the same pattern, so it is enough to actually prove that the re- 
lations (C-4) and (C-7) are true. This can be done in steps as follows: 
a. By trigonometric expansion of the relation 


% 


for: HR ο ο i 


И 


соз (КВ - [Ty sind 1) 


(cosVO cos8 + sinV@ sind) cos[ Tw sind J + 


Ages 


(sinV9 cost - sind со5 Уд) sinl Tw sin] 


Xy cos VO cos Tw sind J + sin sinL Tu sind ] 


=(cos(V+1)8 cos8 + sin(V4+1)0 sind)cosL Tv sind] 4 


(sin(V+1)8 cos8 - cos (V+1)0 sin0)sinl Ty sinó] 
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е 
| 


νε] 7 cos (V+1)9 cos[ Tu sin9] + sin(V+1)9 sinl T% sin] 


l) 


(соѕЅ\№ cosð - ѕіпуӨ sinü)cosl To sin] + 


(sinV соз0 -+ cosv9 sim )sin| Tu sine] 


| 


X042 (cos (V+1)0 cosl Tw sin] 4 sin(V41)8 sinl Tw sind })cos® j 
(cos (V+1)8 зіп ТО sind] - sin(v+1)0 cos[ Tx sind ])sind 

b. By observing that: 

Xl = Х//сов9 - sind sinlW - Tu sind], and 

Ku 
which added together, they give 

X = 
X + ТЕ) 2 со50 Xy er 
τ uo 
= 0 = ° е B ° 

Хо? X 105 sind sin! (V+1)9 Tu sind] and 
Xy, 15050 + sind sinl (V+1)0 - To sind | 


X cose + sind sin W - Tu sind] 


Ху 


[| 


which added together give 


Х |? + Xy = 2 созо νι) OE 
BENE per 
Cx Now since X = -2 X; š ко holds for K = 1 as well as 


A SY, and Y + l, it holds for every other value of K. Q.E.D. 


S The proof that equations (I-12) are also valid follows the 
Same procedure discussed from paragraph 1 up to 4 of this Appendix. But 
in order to prevent unnecessary duplication one can show only some portion 
of the last part of this ανα. for demonstration: 
а. Since T, = cos КӨ апа 0 = ѕіп КӨ Ьу trigonometric ex- 


pansion follows: 


(1) Tey = cosK9Ə cos9 + sinK0 sind 
Ty = coskK9 
T = соз КО со50 - sinK9 sin9 
K+1 
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(2) D = sink cos9 - cosk9 sino 
Ug ZFS inko 
= 1 6 1 Q 
Uk41 sinK9 cosó + cosKB sin 
b. By considering the expressions for To 1 а? апа Un; опе 


ғаш сес: 


Ө + 0 sin 
К-1 T,cos <р Y, Sin Tor 


ES 
i 





-1 
Uk = (sind ) (-соѕ Т, + 1 1), or 
2 25 
ШЫ ae) er m) (с-10) 
С; One can observe also that: 
sind, 
11 = cos9 ( cosK9 + cogo Sin КӘ) 
Bug 
Tea = соѕ0 ( cosK9 - cogo Sin КӘ) 


Therefore, by adding these two equations, the recurrence relation needed, 





follows: 
Т + Te-1 = 2 cos K = οτι, or 
Три = 20T, Е Т1 (CTO 
U._1 = cosd ( Ue . £m соз КО) 
Вон — eos В cert сов КӘ) 
and again by addition it follows: 
Ur = -260, = ET (0-12) 
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APPENDIX E 
SOLUTION TO A QUARTIC EQUATION 


l. A quartic equation 


x + ES + bx” + cX + d (E-1) 
has the resolvent cubic equation 
y? - БУ + (ac-4d)Y - Pd + 4bd - ша = (E-2) 


ДА By using the attached solution to a cubic equation, one can specify 


arbitrarily a solution y = У, to equation (E-2) and then define 
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ATTACHMENT TO APPENDIX E 


SOLUTION TO A CUBIC EQUATION 


ME: Monsider the Cubic: 
3 2 
ΡΥ + qY¥ + r = 0 eo) 
Let: y = x - E C2) 
Calculate: 
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3 27 


and substitute into (E-1) , it follows: 
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